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ABSTRACT 

We study the case of DM self annihilation, in order to asses its importance as an energy 
injection mechanism, to the IGM in general, and to the medium within particular DM 
haloes. We explain the fact that for thermal relic WIMPS, the annihilation cross 
section is model independent (except for possible low velocity enhancements) , and we 
consider two well motivated WIMP candidates, the SUSY neutralino and the first KK 
excited state of the B electroweak boson (also called the KK-photon) , that have masses 
in the range favoured by both, indirect searches and accelerator constraints (m^ > 
lOOGeV). We give an overview of the physical scenarios in which these particles arise, 
mentioning the motivations and general aspects of the underlying physical theories. 
We then analyse in detail the clustering properties of CDM, on all hierarchy levels, 
from haloes to subhaloes and DM density profiles within them. In this analysis, we use 
the Press-Schechter formalism to compute the DM halo mass function, we consider the 
subhalo mass function calibrated from N-body simulations, and we also consider the 
process of adiabatical compression by the presence of a SMBH for our computations. 
We then explicitly compute the energy output (or luminosity) of DM haloes due to 
annihilations, and compare the obtained luminosities with the standard AGN feedback 
process, concluding that DM annihilation does not provide the necessary output as to 
constitute an important feedback mechanism. We then compute the energy injection 
rate per baryon of annihilations on the IGM, in order to calculate the effects that it 
has on its temperature and ionization fraction. We find significant deviations in the 
evolutions of the temperature and ionization fraction of the IGM, in scenarios that 
take into account the clustering of DM at all levels, such that a ITeV WIMP may, for 
example, maintain the temperature of the IGM on the ~ lOOK level, or contribute to 
the ionization fraction maintaining it in a value of ~ 10~^ (not considering POPIII 
stars). We conclude that although the DM annihilation can have a significant impact 
on the properties of the IGM, it cannot be regarded as an alternative reionization 
scenario. Regarding the detectability of the WIMP through the modifications to the 
21 cm differential brightness temperature signal (5Tb), we conclude that, considering 
the adiabatic compression and the substructure, and assuming a maximal absorption 
fraction of fabs = 1, a ITeV WIMP with thermal relic cross section may be detected 
with a peak signal of ^ 18mK at a redshift of ^ 50, corresponding to a frequency of 
30MHz. 
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1 INTRODUCTION 

The Dark Matter (DM) is a component of the total 
density of the universe which can be measured only through 
gravitational interactions and (likely) weak interactions. It 
behaves as dust-like matter in the sense that it is collision- 
less and that its energy density decreases as the third power 
of the scale factor. There are many compelling evidences of 
its existence, although all of them are indirect. These are 



mainly dynamical evidences, which rely on the gravitational 
influence of DM. Some examples of these evidences are the 
rotation curves of late type galax ies (jRubin fc Ford||1970 l, 
galaxy cluster dynamics ( |Zwicky||1933[ ) and strong lensing 
( Paczyhski|1987| ), among others. The elTects of DM are very 
important not only for determining the dynamics of galaxies 
or galaxy clusters, but also for the evolution of the universe 
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as a whole, and for the galaxy formation processes and struc- 
ture formation scenarios. 

For example, assuming that the DM consists of par- 
ticles that were originally in thermal equilibrium with the 
radiation in the early universe, depending on the mass of 
the particles, the structure formation scenario will be top 
down or hierarchical. In the top down structure formation 
scenario, the first structures to collapse are large (on the 
mass scale that corresponds to present-day superclusters), 
because all overdensities at smaller scales are erased by free 
streaming. This is the case for a hot dark matter (HDM) 
model, in which the DM particles were relativistic at the 
time of decoupling. In the hierarchical structure formation 
scenario, small (and also large) structures start to collapse 
after decoupling, because DM particles were non-relativistic 
at the time of decoupling, and so the free streaming is negli- 
gible. This is the case for a cold dark matter (CDM) model 
( Blumenthal||1984 i. 

The postulated CDM particle has been named a WIMP 
(or Weakly Interacting Massive Particle), because it weakly 
couples to normal matter (for introductions to the topic of 
DM astrophysics, se e|D'Amico, Kamionkowski fc Sigurdson] 
1 2009[ ) and |Raffelt| ( |2008[ )). In the standard scenario, this 
particle is considered to be a thermal relic. This means that 
its creation and annihilation processes from the early uni- 
verse radiation field were in thermal equilibrium with each 
other, and thus, the present day DM number density is the 
relic abundance after the freeze-out of this processes, with 
annihilation being favoured after the universe had cooled 
enough ( Kolb fc Turner] 1990 1 . This fact alone allows for the 
computation of the interaction cross section for the creation- 
annihilation process given only the present day abundance 
of DM (its density parameter Q,m,o)- 

There are several DM candidates and scenarios that al- 
low for the existence of a thermal relic WIMP. One of the 
most highly motivated is the Supersymetric Lightest Neu- 
tralino (or simply neutralino), corresponding to the light- 
est mass eigenstate of the neutral supersymmetric part- 
ners of the Standard Model electroweak gauge bosons and 
Higgses (see for exampl e |Martin| ( [2008| ); and IJungman^ 
Kamionkowski & Griest (1996l). In Supersymmetric theo- 
ries (called SUSY in general), there is one boson for every 
fermion and the corresponding partners are called super- 
partners. As the gauge bosons are vector bosons and the 
higgses are scalar bosons, the corresponding superpartners 
are fermions. In SUSY, the neutralino is a Majorana fermion, 
meaning that it is its own anti-particle, and thus it can 
annihilate with itself. Also, in R-parity conserving models, 
the neutralino is stable because no supersymmetric parti- 
cle is allowed to decay into standard model particles. Thus, 
the SUSY neutralino may provide a suitable thermal relic 
WIMP candidate. 

Another highly motivated DM candidate is the Light- 
est Kaluza Klein (LKK) excited state. This corresponds to 
excitations of standard model particles but in a higher di- 
mensional space coordinate. Spatial extra dimensions is a 
prediction of several beyond-standard-model physical theo- 
ries, for example Superstring Theory. In this theories, the 
extra space dimensions are compactified, so that the effec- 
tive theory in the low energy limit is four dimensional. Be- 
cause of the compactified nature of these dimensions, the 
momentum in them is discreet, and in the effective theory 



this extra momentum will take the form of an additional 
contribution to the rest-mass of the particle. KK particles 
are thus heavy, and conservation of momentum in the ex- 
tra dimension guarantees that they are stable. For giving an 
appropriate WIMP, charge neutrality is also required, and 
thus, the best KK candidate for the WIMP is the S'^' (the 
first excited state of the B gauge boson, roughly equivalent 
to the first excited state of the photon). The KK WIMP 
is a boson, whereas the SUSY WIMP is a fermion, however 
both are their own anti-particle and thus can annihilate with 
themselves ( Hooper fc Profumo|2007 |. 

An important feature of a self-annihilating thermal relic 
WIMP is that the annihilation products are standard model 
particles, and thus, they may inject energy to the IGM or 
to the medium where a particular DM halo is embedded. 
Globally, the annihilations of WIMPS are unlikely to oc- 
cur because their abundance has long ago frozen. However, 
there may be regions where encounters among WIMPS are 
more probable due to an enhanced density. These regions 
are precisely the regions where structures form. It is a di- 
rect consequence of the ACDM model that DM is clumpy, 
due to the hierarchical growth of structure. Thus, it is pre- 
cisely in the inner regions of DM haloes and substructures 
that the annihilations are more probable. The energy out- 
put due to annihilations will be proportional to the square 
of the density, thus, these local density enhancements are of 



great importance ( Natarajan, Croton fc Bcrtone^^2008) ). 

The first hierarchy in density enhancements are the 
DM haloes that do not belong to any gravitationally bound 
structure at a particular redshift. In N-body DM only simu- 
lations, these are called main haloes. An analytic formalism 
to quantify the number density of these haloes per decade 



in mass is given by the Press-Schechter formalism (Press fc 



Schechter|(|1974[);|Mo fc White] (|1996|); [Bond et al | ( |1991| ); 
jsheth. Mo fc Tormen| ( |2001[ ); [Lacey fc Cole| ( |1993f among 
others). This formalism considers that the overdensity field 
of the universe is a gaussian field, and thus, the proper- 
ties of its fourier transform are known. One can define the 
matter power-spectrum as the square of the complex mod- 
ulus of the overdensity in k space, and compute the vari- 
ance of the matter fluctuations of scales of a given mass, as 
a window-smoothed average of the power spectrum on the 
corresponding volume. The Press-Schechter formalism then 
allows to obtain the number density of haloes considering 
that the gravitational collapse is a self-similar process, and 
thus, a given halo collapses if its mass scale variance sur- 
passes a certain theshold for collapse at that redshift, and 
the probability that a halo has the required mass is dis- 
tributed gaussianly with the variance corresponding to the 
associated mass. 

The second hierarchy in density enhancements are the 
substructures embedded in a particular halo. As the gravita- 
tional collapse is a self-similar process, each halo will be the 
end node (the trunk) of an inverted halo merger tree at that 
redshift ( [Springel et al.|2005[ ). As the DM is non-collisional, 
it relaxes much slower than normal matter, because it can- 
not radiate away energy to fall to the bottom of the grav- 
itational potential well. Also, as friction is negligible, the 
angular momentum barrier will make the DM subhaloes to 
acquire stable orbits around the center of mass, but these 
subhaloes will retain their identities. Thus, inside each DM 
halo, it is expected that the full hierarchy of DM haloes 
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is nested within, until the lowest mass haloes that formed. 
The number density of these subhaloes can be measured 
in Milky- Way sized halo simulations, and by approximating 
their abundance by a power-law behaviour (Via Lactea II 
simulation (Diemand et al. 20081, Acquarius project simu- 
lation i jSpringel et al.||2008| )) 



example, Cline, Vincent & Xue (20101 and Finkbeiner et 
ah] ( |2010[ ) analyse the data sets by PAMELA, Fermi and 
and conclude that the observed signals could be 
ITeV WIMP. They also find an en- 



Apart from the number density of haloes and subhaloes, 
it is important to quantify the form of the density distribu- 
tion within a particular structure, in order to compute the 
annihilation power output. [Navarro, Frenk fc White| ( |1997[ ), 
when studying the clustering of DM in N-body simulations, 
found that haloes from different Cold Dark Matter cosmolo- 
gies all have the same universal density profile (the NFW 
profile). This profile is characterized by two parameters, the 
scale radius and the scale density, and this parameters can 
be (on average) regarded as a function of the halo mass 
and the redshift. Also, this DM density profile may be fur- 
ther enhanced by the gravitational effect of baryons on DM 



H.E.S.S 

accounted for by a 
hanced interaction cross-section (see Sommerfeld enhance- 
ment. Section 4.5) for the annihilation process (~ lO'^ times 
the thermal relic WIMP value). A similar result was also 



obtained by Elsasser & Mannheim ( 2004 1 with respect to 
the EGRET gamma ray data ( Sreekumar et al.|1998j ). How- 
ever, there is some debate regarding the compatibility of the 
presence of Sommerfeld enhancement with the non-detection 
of gamma ray flux from known DM structures (like galaxy 



clusters). In this line, the study by Pinzke, Pfrommer & 



through adiabatic compression (see, for example Sellwood & 



McGaugh ( 2005 I; and Natarajan, Croton fc Bertone ( 2008 1) 



Bergstrom ( 2009 1 determines that the attempts of explain- 
ing the Fermi, PAMELA and H.E.S.S. combined data with 
~ ITeV WIMPS and cross sections enhanced by factors 
of 10^ - 10^ are incompatible with the EGRET upper 
limits on the gamma ray emission from the Virgo galaxy 
cluster. Regarding the Fermi gamma ray excess, [Hooper 
Goodenough (20101 argue that it can be well accounted for 



by non-exotic processes (decaying pions, inverse Gompton, 
point sources, etc.), except within 1.25° (~ 175pc) from 
the galactic center, where the gamma ray excess is consis- 
tent with the prediction of a cusped (p oc r"^'^'') DM den- 
sity profile, and annihilating WIMPS with a thermal relic 
cross-section (no Sommerfeld enhancement) and masses in 
paoli| (|2006[); and [Gumberbatch, Lattanzi fc Silk] (|2010[)). the range = 7.3 - 9.2GeV. Also, [Hooper, Finkbeiner] 



Thus, the total resulting DM distribution in the universe can 
be computed at a particular redshift. 

Knowing the global DM distribution, one can compute 
the resulting energy output per baryon, which can be further 
compared with what can be obtained assuming a perfectly 



homogeneous DM distribution (see Mapelli, Ferrara & Pier- 



One can also compute the energy generated per unit time 
by a particular halo of a given mass and redshift, and one 
can calculate what is the energy received by a particular 
halo due to the rest of the universe. 

This additional energy injection may have interesting 
effects on the global IGM, as well as on the environment 
within particular DM haloes. For example, several authors 



(for example Gumberbatch, Lattanzi fc Silk, (2010); Mapelli, 


Ferrara & Pierpaoli 


(20061; Ripamonti, Mapelli & Ferrara 


(|2007|);|Valdes et al. 


(|2007|); |Pierpaoh| (|2004|) and|Natarajan 


fc Schwarz (2009l) have studied the effects of DM annihila- 



tion (and also decays, in different DM models) on the IGM 
temperature and on the global ionization fraction. This extra 
heating may be detectable through the global 21 cm differ- 
ential brightness temperature, measurable in the MHz fre- 
quency range. This differential brightness temperature arises 
as a result of the coupling of the neutral hydrogen spin tem- 
perature (that accounts for the hyperfine transition) to the 
IGM kinetic temperature, thus causing a difference between 
the spin temperature and the GMB temperature. Also, the 
energy input may contribute as an extra heating mechanism 
within haloes that can add up to the standard SNe and 



AGN feedback mechanisms (Natarajan, Croton & Bertone 



(20081), and may change the temperature of the halo envi- 
ronment. 

Regarding the observational evidence for the energy in- 
jection due to DM annihilations, no definite conclusion has 
been drawn yet. However, different observations of anoma- 
lous enhancements in detected signals from high-energy phe- 
nomena in the galactic halo and in the cosmic background 
suggest its presence. These signals correspond mainly to 



gamma ray detections (Fermi (Abdo et al. 20091, H.E.S.S 
( Aharonianetal.|2008|) and EGRET (iSreekumar et al.|1998 1 



experiments) and to measurements of cosmic ray positron 
fractions (PAMELA ( [Adriani et al.||2009[ ) experiment). For 



fc Dobler ( 2007 1 found that the excess microwave radiation 
within 20° from the galactic center, found by the WMAP ex- 
periment and called the WMAP Haze, can be well explained 
by a cusped DM profile with p cx r~^'^ in the inner kilopar- 
secs, and a WIMP with mass in the lOOGeV to multi-TeV 
range and thermal relic cross section. It can be seen that 
no clear consensus exists neither on the significance of the 
results as positive indirect detections of DM, nor on the im- 
portance of these detections for discriminating among differ- 
ent WIMP scenarios. What can be said nonetheless, is that 
evidence of possible new exotic physics is accumulating, and 
that WIMP DM annihilation represents a well motivated 
option. As for the evidence for Sommerfeld enhancement, it 
is specially controversial and largely model-dependent. How- 
ever, the DM annihilation interpretation of the observational 
evidence supports WIMP masses in the range ~ lOGeV and 
~ ITeV. 

The main goal of this work is to compute this addi- 
tional energy injection, both locally within haloes and glob- 
ally, and to asses the importance of this injection on the 
different processes mentioned above. All the computations 
will be done considering a cosmology consistent with the 
First Year Data Release of WMAP ( [Spergel et al.|2003[ ). In 
particular, we will use h — 0.72, ^m,o = 0.27, Q.a,o = 0.73, 
fbary = = 0.17, Us = 0.99 and as = 0.9. The paper is 

organized as follows. In Section 2, the thermal relic WIMP, 
and the two above mentioned well motivated particular DM 
models will be explained in greater detail. In Section 3, the 
statistics of the DM clustering on all hierarchies will be de- 
veloped and the NFW universal density profile, as well as 
the adiabatic compression mechanism will be explained. In 
Section 4, the DM annihilation energy output, both globally 
(per baryon) and locally within a particular DM halo will 
be computed. Also, the dumpiness factor, defined as the 
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ratio of output between the smoothed and fuUy clumped 
cases, will be calculated. In Section 5, the energy injected 
on the medium of a particular halo, due to the contribu- 
tion of the rest of the universe, will be obtained for different 
cosmological scenarios. In Section 6, the effects of the anni- 
hilation energy injection on the global IGM will be followed 
in detail, computing the IGM temperature, the ionization 
fraction and the resulting 21 cm differential brightness tem- 
perature in this scenario. In Section 7, the formalism for cal- 
culating the effects of the DM annihilation energy injection 
on individual halos, in particular the change in temperature 
and ionization fraction of the intra halo medium (IGM for 
cluster-sized haloes), will be outlined. The general conclu- 
sions will be presented in Section 8. 



2 WIMP MODELS AND THE THERMAL 
RELIC 

As was previously mentioned in the introduction, the 
WIMP is a theorized particle that allows to account for the 
missing mass in the universe. The evidence of this miss- 
ing mass is (as previously mentioned) mainly gravitational. 
This missing mass cannot come however from normal bary- 
onic matter, because current Big-Bang nucleosynthesis con- 
straints, like for example the primordial Helium abundance 
(Section 10.4 of [Ryden] ( |2003[ )), set stringent limits on the 
normal matter content. 

For a certain particle, motivated by some beyond- 
standard physical scenario, to be a viable WIMP candidate, 
it must have some basic properties. For example, it has to 
be charge-neutral and it can interact only weakly with ordi- 
nary matter. Also, it has to be coUisionless and massive in 
order to behave as Cold Dark Matter. A light DM particle 
will behave as Warm or Hot dark matter depending on if it 
was relativistic or not at the time of decoupling and also de- 



pending on its free streaming length ( Angulo & White 2010 1, 



and the resulting structure formation scenario or minimum 
halo mass will be different to the CDM case. Also, a WIMP 
must be stable, or have a decay time longer than (or at least 
comparable to) the present Hubble time. 

In some well motivated physical scenarios, the WIMP is 
a thermal relic. This means, as was mentioned before, that 
its creation and annihilation reaction were in thermal equi- 
librium with the radiation background in the early universe. 
Then, as the universe cooled down, the annihilation reaction 
was increasingly favoured over the creation reaction due to 
the rest-mass energy difference between both states. Finally, 
as the expansion rate of the universe became larger than the 
annihilation reaction rate, the global abundance of DM (its 
comoving number density) was frozen, and the probability 
of occurrence of an annihilation event became negligible. 

In two well studied scenarios, the WIMP particles are 
such that they are their own anti-particle. For example in 
the case of a Kaluza-Klein excited state, arising naturally in 
compactified extra-dimensional models (like Universal Ex- 
tra Dimensions or UED, in which all particles are able to 
propagate in the extra dimension provided they have enough 
energy), the WIMP is the lightest neutral gauge boson ex- 
cited state, corresponding to the first excited state of the 
B boson or B^^\ and as the standard model B boson, it is 
its own antiparticle (the B boson is a superposition of the 



photon of Quantum Electrodynamics and the Z boson of the 
weak interaction). Also, in the case of the MSSM (the min- 
imal supersymetric standard model), considering conserved 
R-parity (and thus that superpartners are forbidden to de- 
cay to standard model particles), the WIMP is the lightest 
neutralino (i.e., the lightest mass eigenstate of the fermionic 
superpartners of the gauge bosons and the Higgses), and it 
is a Majorana fermion, and thus it is its own antiparticle. 

In both these scenarios, the annihilation cross section is 
dependent only on the present day abundance of DM. This 
fact has been dubbed the WIMP miracle and is explained in 
what follows. Also, both scenarios are explained in greater 
detail. 



2.1 The WIMP miracle 

The WIMP miracle refers to the fact that, given that 
the WIMP is a thermal relic, its creation/annihilation cross 
section can be computed without needing extra model- 
dependent input (except for the fact that the cross section 
must not depend on the energy of the WIMP, i.e., the s- 
wave cross section must be dominant). The cross section 
then only depends on the present day DM density parame- 
ter Qm,o, and the present day value of the Hubble parameter 
h (that gives the current expansion rate). 

Intuitively, one would think that a higher cross section 
< av > would lead to less relic abundance of DM, thus 
resulting in a smaller value of Qrn.o- Both quantities are, 
indeed, inversely proportional to each other. 

Obtaining this cross section is a standard derivation 
and it is done in, for example. Section 5 .2 of|Kolb fc Turner 
(|1990|). For the value of the cross section. Hooper, Finkbeiner 
(2007[) use < cru >= 3 X IQ-^^cm-' / s], and in the 



& Dobler 



cosmology used in this work, that gives: 



< av >- 



4 X 10"^ 



(1) 



In all that follows, we will use that value for the cross 
section of WIMP annihilation. 

Now we will motivate and explain two scenarios that 
give rise naturally to a thermal relic WIMP with the above- 
given cross section. 

2.2 The KK DM particle 

In extra dimensional models (i.e., physical models that 
are built considering more than three spatial dimensions), 
the additional spatial degrees of freedom may be accessible 
to all or some of the particles. Also, these additional dimen- 
sions may be compactified differently and on different energy 
scales. If a standard model particle can propagate in this 
compactified (and thus finitely extended) dimensions, they 
will appear more massive in the four dimensional effective 
theory. In some cases, these excited states of standard model 
particles may have the desirable features of the WIMP. Here 
we will briefly discuss the lightest KK particle WIMP can- 
didate and the extra-dimensional scenario in which it arises. 
We loosely follow the excellent review by |Hooper fc Profumo| 
( 2007] ). 

The first one to seriously propose that there may be 



more spatial dimensions was Theodor Kaluza (Kaluza 1921 1, 
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in an attempt to explain electromagnetism in the framework 
of General Relativity. Although this proposal was not con- 
sidered by mainstream physics, recently it was demonstrated 
that string theory could not be formulated in only four di- 
mensions, and thus, this idea was again taken into account. 

One extra dimensional scenario that is interesting for 
being able to provide a suitable DM particle is Universal 



Extra Dimensions (UED). In this scenario, proposed by Ap- 



pelquist et al. (2001 1, all standard model (SM) particles are 
able to propagate in the extra dimension. Because this di- 
mension is compactified, and thus has a finite extent, the 
momentum that a particle can acquire in it will be discreet 
(the idea is similar to that of an electron propagating in a 
box of given size), thus giving an apparent extra rest-mass 
in the four dimensional world. 

The spectra of different possible momentum states in 
the extra dimension will give a series of possible masses of 
a particle, quantized by the number of the excitation state. 
This masses will depend on the radius of compactification 
and will be given roughly by: 



i?2 



-I- m 



(2) 



where m_y(„) is the mass of the n' excited Kaluza-Klein 
state of a SM particle X, m^{o) is its mass in the SM and R 
is the compactification radius of the extra dimension, given 
in natural units (so it has dimension of Energy'^). Then, 
the KK particles are the same SM particles, but consid- 
ering their excitations in the extra dimension. Interactions 
between SM particles and KK particles are possible, but 
are suppressed (thus, these particles are weakly interact- 
ing). Moreover, the lightest KK particle is stable because of 
conservation of momentum in the extra dimension. 

As was already mentioned, the lightest KK particle in 
the UED framework is the S'^', i.e., the first KK excited 
state of the B gauge boson (also called the KK- photon). This 
particle is charge neutral, and can annihilate with itself, and 
thus is a suitable WIMP candidate. 

Some of the ways in which KK particles can interact 
with SM particles are through decays, annihilations and 
scatterings. For example, a heavier KK particle can decay 
into a lighter KK particle emitting standard model particles, 
until it has cascaded into the S'-^-' (also called 71, see Figure 
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Also, more interestingly in view of the possible effects of 
WIMP annihilation on the IGM, two S'^'s can annihilate 
to lepton pairs and photons (see Figure 16 of [Hooper fc 
|Profumo| ( [20071 )). 

Finally, the KK-particles and in particular, the S'^-*, 
can scatter off SM quarks and leptons (see Figures 10 and 
12 of [Hooper fc Profumo[ ( [2007| ). Ahhough, these processes 
are highly suppressed (making the S'^' a weakly interacting 
particle), they can occur, and so, this WIMP candidate can 
be searched for using direct detection methods (for example, 
scatterings off heavy nuclei). 

Although of importance for particle physics and for con- 
straining the UED parameters in collision experiments, a 
computation of the branching ratios and cross sections of 
the processes is beyond the scope of this work, and it is also 
not required for the computations performed in this work, 
because the resulting WIMP is a thermal relic. 



2.3 The Neutralino 

For this section, we will follow two reviews that analyse 
different theoretical aspects of the SUSY neutralino. The 
first review, by iMartin' (2008), gives a thorough introduc- 
tion to the MSSM, defining what is supersymmetry, giv- 
ing the MSSM lagrangian, explaining the different SUSY- 
breaking scenarios (that ultimately fix the parameters of 
the particular model), etc... The second review, by Jung- 
[man, Kami onkows ki fc Griest[ ( [1996[ ) , treats the neutralino 
DM WIMP candidate, and its associated physical features, 
with in-depth detail, calculating the particular cross sec- 
tions of all the interaction channels of the neutralino with 
SM particles for the most common SUSY scenarios. 

Supersymmetry is a particular symmetry property of 
the lagrangian of a physical theory, such that said lagrangian 
remains invariant if a supersymetric transformation is per- 
formed on the fields of the theory. The supersymmetric 
transformation is such that the change in a fermionic field 
is a bosonic field, and the change in a bosonic field is a 
fermionic filed. As in quantum field theory, the fields rep- 
resent different particle contents (the particles are nothing 
more than the excited states of the fields), a supersymmet- 
ric theory naturally has one boson for each fermion present, 
and vice-versa. 

The lagrangian of the SM of particle physics is not su- 
persymmetric. As the SM has been remarkably successful 
in accounting for all known physical interactions (except 
for gravity), a natural question is why would it be needed 
to propose a supersymmetric theory. As of today, the two 
main motivations for invoking supersymmetry (SUSY) are 
the unification of the interaction coupling constants at high 
energy and the solution of the hierarchy problem. 

It is a known fact in quantum field theory, that the cou- 
pling constants are a function of the energy (or momentum) 
of the interaction. This is called the running of the couplings 
and it is a consequence of the infinite loop contributions to 
any given process. Although a more detailed explanation of 
this phenomena is beyond the scope of this work, it is suf- 
ficient to say that by using renormalisation techniques, it is 
possible to find the dependence of the couplings with the en- 
ergy, and that when considering SUSY, the couplings tend 
to unify at an energy scale of IQ^^GeV (see Figure 5.8 of 
Martin[([2008|)). 



The strong hierarchy problem comes from the quantum 
corrections to the lagrangian mass of the Higgs field in the 
standard model. What happens is that for every particle in 
the model, the Higgs mass receives a correction that is lin- 
ear on the energy scale at which the SM ceases to be valid. 
Of course, if the SM is the ultimate physical theory, then 
this would not be a problem, however, this cannot be be- 
cause the SM does not explain gravity. Reasonable order of 
magnitude arguments point in the direction that, at least 
at an energy scale of the order of the Planck Mass (Mp), 
quantum-gravitational effects should be important (for ex- 
ample, the fact that for a particle of mass Mp, its Compton 
wavelength is comparable to its Schwarzschild radius), and 
so, the SM ceases to be valid at least at this energy scale. 
SUSY naturally solves this problem, because the contribu- 
tion to the Higgs mass is positive in the case of bosons and 
negative in the case of fermions, and so as for each boson 
there is a fermion, the contribution cancels out exactly. 
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The particle content of the MSSM (the Minimal Super- 
symmetric Standard Model) is the same as that of the SM, 
but for every SM particle, there is an associated supersym- 
metric partner (also, the Higgs field is different from the one 
in the SM, being two fields instead of just one). The naming 
scheme for the superpartners is to maintain the root of the 
name of the SM particle, but to change the end of the name 
with the suffix -ino in the case of fermionic superpartners of 
SM bosons and to add the prefix s- in the case of bosonic 
(scalar) superpartners of SM fermions. For example, the su- 
perpartners of the electron and the muon are the selectron 
and the smuon, and the superpartners of the photon and the 
Z boson are the photino and the Zino. 

As it was mentioned above, the WIMP candidate of 
SUSY is the lightest neutralino (unless the gravitino, the 
superpartner of the graviton in supergravity is lighter than 
the neutralino, but this scenario will not be considered here). 
The neutralinos are the mass eigenstates of the superpart- 
ners of the neutral gauge bosons and the Higgses, i.e., the 
Bino (B), neutral Wino (W°) and the neutral Higgsinos ( 
H° and H^.)- They are labelled Ni,i = 1...4 in order of in- 
creasing mass, such that the A''i corresponds to the lightest 
neutralino. The A'^i, (the WIMP candidate) will be in gen- 
eral a linear combination of the gauginos (Bino and Wino) 
and the Higgsinos, and it will have a somewhat different 
phenomenology depending on its specific particle mixture. 
Although the precise differences between a mainly gaug- 
ino, mainly Higgsino or mixed gaugino-Higgsino type neu- 
tralino are beyond the scope of this work, it is sufficient 
to mention that, in general (in order to give the appro- 
priate above-mentioned interaction cross section), a mixed 
gaugino-Higgsino composition is assumed. 

In R-parity conserving SUSY scenarios, no superpartner 
can decay to SM particles only (SM particles have different 
R-parity than superpartners) , and so the lightest neutralino 
(A^i) will be stable, whereas more massive SUSY particles 
can decay into the Ni emitting standard model particles (as 
long as the standard conservation rules are not violated). 
Also, the (A^i) is a Majorana fermion, and thus it can self- 
annihilate into SM particles. The SUSY particles can also in- 
teract (scatter off) standard model particles, but these scat- 
terings are suppressed, thus the Ni is weakly interacting. 
The main annihilation processes of neutralino into SM par- 
ticles is into fermion pairs, although it can also annihilate 
into gauge bosons (W and Z) at tree-level. At 1-loop level, 
the neutralino can annihilate directly into gluons or pho- 
tons, and although this is greatly suppressed compared to 
tree-level processes, the direct decay into photons is impor- 
tant because it signals a smoking gun of annihilation (the 
resulting two photons will have an energy very close to the 
neutralino rest-mass energy, and so it is a clear signal with 
no contaminants). The Feynman diagrams of the tree-level 
annihilation channels are shown in Figures 9, 10 and 16 of 



Jungman, Kamionkowski & Griest ( 1996 1 



Like in the case of the lightest KK particle, the Feyn- 
mann diagrams were included just to illustrate the available 
annihilation channels, and a computation of the precise cross 
sections for each process is beyond the scope of this work (al- 
though it is done explicitely by [Jungman, Kamionkowski fc| 
Griest I ( p96t ). 

Finally, it is important to mention that the phe- 
nomenologies of KK particles (as they arise in UED) and of 



SUSY are very similar, and both of them could be tested, for 
example in particle colliders, by studying interactions above 
a certain energy scale. In the case of KK, this energy scale 
corresponds to the inverse of the radius of compactification 
of the extra dimmension, whereas in SUSY, it corresponds 
to the mass of the neutralinos and superpartners. Current 
collider lower limits on the mass scale of this new physics 
(for example the LHC limits) are on the order of the 100 
GeV and increasing. This mass limit will motivate the val- 
ues that we adopt for the WIMP mass in the rest of this 
work. We use both a WIMP mass of 10 GeV and of ITeV. 



3 THE DM CLUSTERING 

In standard FRW (Friedmann-Robertson- Walker) cos- 
mology, the universe is assumed to be homogeneous and 
isotropic on large scales. However, the matter is distributed 
in a random but not homogeneous manner on small scales, 
and thus the universe exhibits structure (like voids, clus- 
ters of galaxies, cosmic filaments, etc.). The hierarchical 
structure formation scenario is a natural consequence of the 
KCDM cosmology, and it predicts that the DM accretion 
is a clumped process. In this scenario, structures of all mass 
scales are allowed to collapse, and so the densest and small- 
est objects collapse first, corresponding to small peaks on 
top of larger peaks of the matter overdensity field. In or- 
der to account for the formation of large-scale structure, 
the DM needs to be followed in the gravitational collapse 
process, but considering the corresponding effect of A (the 
cosmological constant) on the background expansion. 

There are two main approaches for following the grav- 
itational collapse of the DM. The first approach is to sub- 
stitute the continuous overdensity field by DM point masses 
and to perform computational N-body simulations of their 
mutual interactions. The second approach is to make use 
of the statistical properties of the overdensity field, and de- 
velop a formalism that allows to analytically compute the 
average clustering of the DM at different redshifts. 

As it was mentioned in the Introduction, the gravita- 
tional collapse process is hierarchical and self-similar, and so 
different levels of clustering are expected. In the following 
subsections we explain the prescriptions adopted to account 
for the DM clustering at all these levels, considering also 
the distribution of DM within a particular halo or subhalo. 
In Section 3.1, we follow the analytic approach, using the 



Press-Schechter formalism (Press & Schechter (19741; Mo 



fc White] (p96|); [Bond et ai.| ( |1991| ); [Sheth, Mo fc Tormen 



(20011; Lacey & Cole (19931 among others) to compute the 



number density of DM haloes per decade in mass. In Section 
3.2, we account for the presence of substructure by consider- 
ing the substructure mass function given by |Giocoli, Tormen| 
& van den Bosch (20081. In Section 3.3, we present the uni- 
versal NFW density profile ( [Navarro, Frenk fc White|1997'l ), 
that gives the DM distribution within each gravitationally 
bound halo, and explain how to calculate its parameters for 
different halo masses at different redshifts. Finally, in Sec- 
tion 3.4, we explain the effect that the gravitational collapse 
of the barionic matter (in particular, the formation of a su- 
permassive black hole (SMBH)) can have on the DM density 
profile through adiabatic compression. 
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3.1 The halo mass function 

The analytic approach to obtaining the clustering 
statistics of DM has been developed in time by different 
authors. It was in |Press fc Schechter| ( |T974| where a formal- 
ism for computing the halo mass function was first obtained, 
starting from the power spectrum of matter density fluctu- 
ations. Later on, this formalism was refined by |Bond et al.| 
(19911, where it was cast in a more rigorous mathematical 



form (excursion set theory and random walks). Also, it was 
extended by [Lacey fc Cole| ( |1993[ ) for computing the halo 
merger rates at different redshifts, and also the average num- 
ber of mergers a halo has undergone since its formation as a 
function of mass and redshift (they also give a prescription 
to generate montecarlo merger trees). A further extension of 
the formalism was done by Mo & White ( 1996 1, to compute 



the clustering statistics (correlation functions and bias) of 



haloes. Finally, Sheth, Mo fc Tormen (20011 slightly modify 



this formalism (in particular, the functional form of the halo 
mass function) to account for ellipsoidal collapse due to the 
non-sphericity of overdensities (for a discussion of triaxial- 
ity of DM haloes, see |Piffl fc Sotnikova| (|2008|)). We use this 
formalism, as presented in Mo fc White (19961, but using 



the modified Sheth, Mo fc Tormen ( 2001 1 halo mass func- 



tion and bias factor, calibrated with the GIF simulations 
( Kauffmann et al.|[ll99 l. 

We consider the overdensity 5{x) defined as 5{x) — 
, and the overdensity threshold for collapse 5c{z), de- 
fined as the density contrast (in the linear approximation) 
required for an overdensity in a certain region to already 



have formed a collapsed object (Navarro, Frenk & White 



( 1997 1, equation A14). We consider also the mass fiuctuation 
(j{R), defined as the standard deviation of the matter over- 
density field when smoothed on a scale of size R. Finally, the 
dimensionless mass parameter is defined as = l^xWjMlyj' 
and it is a measurement of the mass of a collapsed region 
(DM halo) relative to the mass of the structures that have 
recently collapsed at redshift z. 



We will use the |Sheth, Mo fc Tormen| ( |2001[ ) probabil- 
ity distribution /{v). The mass fluctuation (j{R) is com- 



puted foUowing iLacey & Colel (119931) and iMo fc White 



(19961, starting from the power spectrum of matter over- 
densities, P{k), at redshift z. When calculating the mass 
function we adopt a primordial power spectrum of the form 
P{k) = Ak"' , where ^ is a normalization. We use the 
WMAPl ([Spergel et ah] (|2003|) scalar index, = 0.99. 



The effects on the power spectrum due to the gravita- 
tional collapse of the DM on sub-horizon scales after the 
radiation-matter equality epoch and the Barion Acoustic 
Oscillations (BAOs) are encoded in the transfer function. 
We adopt the functional form given by Bardeen et al. (the 
BBKS transfer function), which ignores the BAOs. 



3.2 The substructure mass function 

The second level of clustering of the DM corresponds 
to the self-bounded haloes that retain their identity within 
bigger haloes. The presence of these substructures is a natu- 
ral consequence of the hierarchical structure formation pic- 
ture, because haloes at all redshifts are formed by the ag- 
gregation of smaller parent haloes. As the DM is collision- 
less, DM haloes are significantly less disrupted than barionic 



matter during the merging processes, and thus subhaloes 
may survive up to the smallest scales (due to the negligible, 
but nonzero, frcc-strcaming of the WIMPS). Many authors 
(for example |Kamionkowski, Kous hiappas fc Kuhlen (20101; 
Giocoh, Fieri, Tormen & Moreno (2009); Giocoh, Fieri fc 



Tormen] (|2008[); Taylor fc Babul] (|2004j); |Fieri, Bertone fc 



Branchini (2008P) have investigated the clustering proper- 



ties of the subhaloes by means of fitting an analytic form to 
the subhalo mass function obtained in high-resolution nu- 
merical N-body simulations and re-simulations. 

We use the substructure mass function proposed by lGio-] 
coli, Tormen fc van den Bosch (2008), because it also con- 



siders the evolution of the substructure in time, due to the 
combined effects of gravitational heating and tidal stripping 
in the potential well of the main halo. These effects will 
tend to erode the subhaloes, which are the remnants of the 
haloes accreted by the host halo. As was found by these au- 
thors, the mass loss (to the smooth component of the main 
halo) of subhaloes can be approximated by an exponential 
decay of the subhalo mass on a characteristic timescale that 
is proportional to the dynamical time of the main halo. 

Explicitly, the unevolved subhalo mass function (that 
considers all subhaloes with the mass they had when they 
were accreted) is universal, and is given by: 



dN 



dln{mv / Mo) 



Nox- 



-6.283a;2 



aMn 



(3) 



with a = 0.8 and No = 0.21. These values were calibrated 



( 2004 1 simmulations, as well as resimulations done by Dolag 



et al. 



using the GIF Kauffmann et al. ( 1999 1 and GIF2 Gao et al 



(2005) 



Here, m„ represents the unevolved subhalo 
mass (the mass that the halo had when it was accreted) and 
Mo represents the present-day host mass. 

To account for time evolution, the authors find the fol- 
lowing relation between the mass of a subhalo at time t 
(given that it was accreted at time tm), and its unevolved 
mass: 



msb(t) = m^exp 



t tfjl 



(4) 



where msb is the evolved mass and t{z) is the characteristic 
mass loss time, given by 



-{z) = To 



TO 


a{z,Mo) 


-1/2 


\Hiz)] 


-1 


o{z = 0,Mo) 




[ Ho \ 





(5) 



and To — 2.0Gyr. 

In this work, it was useful to approximate the accretion 
time of all haloes by the time at which the main halo col- 
lapsed, corresponding to a redshift Zcoii, which we take as 
the redshift at which it acquired half of its mass (see Section 
3.3). In terms of the look-back time to redshift z {LBT{z)), 
equation 15 can be approximated by: 



rrisi, — rrivexp 



LBT(zcoh)- LBT{z) 
r{z) 



(6) 



Now, the required subhalo mass function is the one that 
includes the mass loss of subhaloes in time, so it should 
be — 4^ j"^" , using the chain rule. Expressing the 
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subhalo mass function in terms of msb instead of m„, one 
obtains; 



With the obtained Zcoii{M, z), we proceed to determine 
po and rs using the following relations: 



dN 
dm 



No 
m 



mK-\z,MH) 



(7) 



where m is the present day (at redshift z), evolved mass of 
a subhalo, Mh is the mass of the host halo and K is defined 
by 



K(z,Mh) = exp 



LBT{z,M{Mh)) - LBT{z) 



r{z) 



(8) 



It can be seen that the behaviour of the subhalo mass 
function is that of a power law with an exponential cutoff, 
and that all mass scales are displaced as a function of time 
since halo formation. The displacement is such that a present 
day subhalo of mass m, corresponded to a halo of greater 
mass (by a factor K), and so its abundance is decreased 
because it is the abundance corresponding to the greater 
mass at the time of halo formation. 

These last two equations give the number of subhaloes 
present in a halo of mass M at redshift z, that have a mass 
between m and m + dm, and so correspond to the subhalo 
mass function that we need. In the next section, we explore 
the smooth density profile of a particular halo, in order to 
account for the last level of the DM clustering hierarchy. 

3.3 The NFW universal density profile 

The last hierarchy corresponds to how the DM is 
smoothly distributed within a particular halo or subhalo. 
This DM distribution appears to be universal, as found by 



Navarro, Frenk & White ( 1997 1 using N-body simulations. 



and is given by the NFW density profile: 



p(r 



(r/r,)(l + r/r.)2' 



(9) 



where po is the characteristic density and is the scale 
radius. Thus, the NFW profile depends on only two param- 
eters, which on average only depend on the mass of a halo 
and the redshift at which it is observed. Before obtaining 
these parameters, it will be necessary to compute the red- 
shift at which a halo that currently (at redshift z) has a mass 
M was formed. 



For this quantity Navarro, Frenk & White ( 1997 1 use 



the redshift at which a halo which currently (at redshift 
z) has a mass M has just acquired half of its mass. The 
authors make use of the gaussianity of the overdensity field 
to compute the redshift of collapse as: 



erfc 



5c{Zcou{M, z)) - Sc{z) 
2[cr2(0.01A'/,2) -cr2(M,z)] 



1/2. (10) 



The interpretation of this expression is that, if we con- 
sider a halo of mass M at redshift z, Zcoii{M,z) will corre- 
spond to the redshift at which a collapsed substructure of 
mass O.OIM was already part of the halo under consideration 
with a probability of 1/2. The above equation is an implicit 
definition of Zcoii{M, z), and the election of the mass of the 
small substructure (of O.OIM) is motivated by comparisons 
to N-body simulations considered by the authors. 



Po = pcrit{z)So{M, z), 



(11) 



where pcrit{z) is the critical density of the universe at red- 
shift z, and So{M, z) is a characteristic overdensity (or den- 
sity contrast) that depends on the mass and redshift of the 
halo whose NFW profile we want to obtain. 
Also, 



\ir{M,z) 

c{M,z) 



(12) 



where r^ir{M, z) is the virial radius of a halo of mass M at 
redshift z, corresponding the the radius at which the average 
density within is equal to the density required for the halo to 
be a self-gravitating virialized structure (in the case of the 
Einstein-De Sitter cosmology, this density is 200 times the 
mean density of the universe). c{M,z) is the concentration 
parameter of the halo. 

The virial radius is given by |Navarro, Frenk fc White] 



(19971 as: 



ryir{M,z) = 1.63 X 10' 



M 



1/3 



m,0 



1 -1/3 



^m{z) 



(1 + z) ^kpc, 



(13) 



where 57m (z) is the matter density parameter at redshift z, 
and the virial radius is given in proper kpc. 
The characteristic density is given by: 



do = 3.41 X 10 \lm{z) \ 



(14) 



Finally, the concentration parameter is implicitly given 
in terms of the characteristic density contrast by the rela- 
tion: 



5o{M,z) 



200 c^{M,z) 
ln{l + c{M,z)) 



c{M,z) 
(l+c(M,z)) 



(15) 



Thus, the density profile of a particular halo or subhalo 
is, on average, completely specified by its mass and redshift. 

In the next section, we proceed to study the conse- 
quences of the collapse of the baryons on the density profile 
of the DM haloes, and obtain a modified density profile for 
the latter considering the growth of a central SMBH, a pro- 
cess that appears to be ubiquitous in galaxy formation. 



3.4 Adiabatic compression 

We consider the modification in the NFW DM density 
profile resulting from the gravitational collapse of baryons. 
Baryons behave dynamically differently from DM in their 
collapse process because they can heat and radiate away 
their energy, as they are not coUisionless and have nonzero 
internal pressure. Thus, as they lose energy, they fall towards 
the center of the gravitational potential well (this process is 
also called baryonic infall or cooling) . This is the main mech- 
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anism that allows the formation of galaxies within the DM 
haloes, and it is also responsible for the adiabatic compres- 
sion of the DM halo density profile. 

Different authors have studied the adiabatic compres- 
sion mechanism, and have given algorithms to compute the 
compressed DM density profile, given the original, uncom- 
pressed profile and the newly grown density profile of the 
collapsed baryonic component. Two of the most used meth- 
ods are the Blumenthal algorithm ( |Blumenthal||1984| and 
the Young algorithm (Young 19801. The Young algorithm 



assumes that the potential well has spherical symmetry, and 
so it corresponds to a central potential problem. The Blu- 
menthal algorithm further assumes that all the orbits of the 



DM particles are circular. Sellwood & McGaugh ( 2005 1 con- 



clude that the Blumenthal algorithm overpredicts the com- 
pression in the inner regions. 

We use a different algorithm and account only for the 
adiabatic compression due to the formation of a point mass 
(i.e., a SMBH) in the center of the potential well. We im- 
plement this algorithm as described in |Natarajan, Croton fc| 
Bertone|(l2008|. 



This algorithm assumes that the density profile is of the 
initial NFW form up to the radius of gravitational infiuence 
of the SMBH (defined as the radius up to which it is the 
main contributor to the contained mass), and that within 
this radius, the density profile is described by a power-law 
with a different (steeper) exponent. It also assumes (like 
the Blumenthal algorithm), that the original orbits of DM 
particles in the NFW profile were circular. 



According to Natarajan, Croton & Bertone (20081, the 



value of the exponent of the density profile in the inner re- 
gion is given by 



spheroid velocity dispersion and BH mass, as presented in 



Murray, Quataert & Thompson ( 2005 I 



Mbh = 1.5 X 10** 



'sph 



200[fcm/s] 



(19) 



The mass of the black hole depends, on average, only 
on the mass of the host halo and on the fact that the host is 
a main halo or a subhalo. The relation between both masses 
IS given by ( [Lagos, Cora fc Padilla||2008[ ): 



2.1 



(20) 



2.9 



(21) 



for the case of a main halo, and 

for the case of a subhalo. 

Now we have all the ingredients for computing the adi- 
abatically compressed DM density profile for a given halo of 
mass M at redshift z. 

Putting all the pieces together, the adiabatically com- 
pressed density profile of a halo has four different regimes, 
and is given by: 



pDM{r) = pNFw{r),r > rsi 



(22) 



PDM(r) = Pnfw(tbh) 



T-BH 



, r plateau < T < TBH, 

(23) 



'Jspike — 2 -|- 



-7 



(16) 



The derivation of the value of the exponent can be found 



Quinlan, Hernquist & Sigurdsson (19951. 



Finally, we need the maximum density and the radius 
of gravitational infiuence of the BH to have the complete 
density profile. For the maximum density, we use: 



max 

Pdm 



tspike < OV > 



(17) 



where is the mass of the WIMP particle, and < ov > 
is the thermal relic annihilation cross-section (see Section 
2.1). For t spike, the time elapsed since the spike was formed, 
we use tspike ~ LBT{zeoii), i.e., the look-back time to the 
redshift at which the halo collapsed (calculated as mentioned 
in Section 3.3). 

For the gravitational influence radius of the black hole. 



TBH = 0.2 



GMbi 



sph 



(18) 



where Mbh is the black hole mass, and a^ph is the velocity 
dispersion of the spheroid component (the bulge in the case 
of haloes typical of late type galaxies). 

For asph, we use the observational relation between 



pDM{r) 



tspike < UV > 



,4rs <r < Tplate 



pDM{r) = 0,r < 4rs 



(24) 



(25) 



where Vpiateau is the radius at which the density of the spike 
equals the maximum attainable DM density, and is the 
Schwarzschild radius of the BH. In Figure 1, we show the adi- 
abatically compressed density profile of a DM halo of mass 
5 X lO^Afg at redshift z = l (Typical of QSO systems). The 
four regimes are clearly distinguishable. The outer slope cor- 
responds to the NFW 7=1 and the inner slope corresponds 

to ^spike — ^ . 



4 DM ANNIHILATION ENERGY OUTPUT 

Having considered the clustering properties of the DM 
at all scales and hierarchy levels, and having explained two 
particular DM WIMP candidates, motivating their expected 
masses and annihilation cross-sections, we can proceed to 
compute the expected energy output of annihilations in lo- 
cal structures (halos) and as an average for the whole IGM. 
In this section, we obtain the luminosity due to annihilations 
of individual DM haloes, considering the clustering analysed 
previously, and we also compute the annihilation energy rate 
injected to the IGM, per baryon, in the universe as a whole. 
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Figure 1. The adiabatically compressed density profile of a halo 



4.1 The smooth DM annihilation energy injection 
rate per baryon 

First we consider the global case, assuming a perfectly 
smooth DM density. This case is treated, for example, by 



Mapclli, Fcrrara & Pierpaoli (20061; Ripamonti, Mapelli & 



Ferrara (20071; and Cumberbatch, Lattanzi & Silk (20101 



As the universe may be infinite in spatial extent (as would 
be the case for flat, i.e. ACDM and Einstein-De Sitter, and 
Open cosmologies), the relevant quantity is the energy in- 
jected to the IGM per baryon. In terms of the present day 
WIMP number density, this is calculated as: 



}.II^c\lo{l + zf <av> Us{z), (29) 
2 Ubfi 



where n^^o is the present day DM number density, and is 
calculated as 



PcO 



(30) 



where ildm.o = (1 ~ fbary)^m,o IS the present day DM den- 
sity parameter, p^.o is the present day critical density of the 
universe and is the WIMP mass. Also, n^fi is the present 
day baryon density, and is calculated as 



We also mention the observational evidence of this annihi- 
lation energy injection process in satellite haloes and in the 
Milky Way (MW). Finally, we asses the importance of the 
DM dumpiness in the energy output, comparing the case of 
a perfectly smooth universe with one with clustered DM. 

First, let us consider the generic energy output. The 
reaction rate of the annihilation process is given by 



riy < av >, 



(26) 



and thus, the variation of the number density of WIMPS 
can be calculated as: 



dn 
'dt 



1 2 

-n^ <av> 



(27) 



As in the annihilation process the WIMPS are de- 
stroyed, we can consider that all their rest-mass is available 
as part of the energy output (actually, only a fraction of 
this energy will be available for affecting the environment, 
depending on the mechanism that couples the annihilation 
products and the IGM or halo medium). Then, the power 
density of the energy output can be written as: 



2 2 

-rriyC < av > 



(28) 



This equation is completely general, and so it applies for 
haloes and also for the diffuse IGM. Let us now consider the 
power output in the different situations (we will call them 
local and global). 



nb.o — ^b,o 



Pc,0 

finiH ' 



(31) 



where ^Itfl is the present day baryonic density parameter 
(equal to fbary^m,o), tuh is the hydrogen mass and /i is 
the mean atomic weight of the baryon content. The mean 
atomic weight, considering a universe with only hydrogen 
and helium (a good approximation, particularly before star 
formation), can be computed as fj. = fn +4/_f/e, where fn is 
the fraction of hydrogen and /hb is the fraction of helium by 
number. Assuming a value of = 0.74 and Y — 0.26 (the 
hydrogen and helium fractions by mass, respectively), gives 
fn = 0.92 and fne = 0.08. Thus, fi = 1.24 will be used in 
this work. 

Knowing the absorption fraction of the energy fabs{z), 
we can readily compute the energy injected to the IGM per 
baryon through DM annihilations. Different authors use var- 
ious prescriptions for /a6s(z). 'Natarajan & Schwarz (20091 
consider that only the photons resulting as secondary annihi- 
lation products can inject energy to the IGM, and calculate 
the photon energy spectrum for different neutralino mod- 
els explicitly. They consider that as the photons propagate, 
they loose energy, and compute the probability that they 
scatter off the IGM atoms. They find that fabs[z) is in the 
range of 0.1 — 0.2 for 2 < 50 (they do not consider higher 
redshifts). 



Cumberbatch, Lattanzi & Silk (20101, compute fabs{z) 



for the neutralino (see also (Ripamonti, Mapelli & Ferrara 



2007 1), calculating first the number of photons and electrons 



produced per neutralino annihilation event, and they obtain 
fabs{z) in the range of 0.01 — 0.1, being higher for higher z 
(and considering redshifts up to z = 1500). 



Other authors, like Mapelli, Ferrara & Pierpaoli (20061, 



simply assume that fabs[z) = 1, and compute the maximal 
effects that their WIMP candidates can have on the IGM. 
We will consider that fabs is a constant that does not de- 
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Figure 2. The smooth energy injection rate per baryon E-^{z), 
for different WIMP masses. The dotted curve corresponds to a 
WIMP mass of ITeV and the dashed curve to a WIMP mass of 
lOGeV. 



pend on redshift, and we will use the values fats ~ 0.01, 0.1 
and 1, in order to obtain results comparable to the different 
authors. 

In Figure 2, we show the results obtained for the smooth 
energy injection per baryon, for two different WIMP masses 
(ITeV and lOGeV), as explained in the caption. We only 
consider the maximal absorption fraction {fabs = 1), be- 
cause for other values, the resulting curves should simply be 
rescaled by fats- It can be seen that the effect of a lower 
WIMP mass is, as expected, to boost the energy injection 



rate by a factor oc 



4.2 The DM annihilation luminosity of a halo 

Having considered the global case of the energy injec- 
tion to the IGM, we now consider the energy output per 
unit time (or luminosity) of a particular halo of mass M at 
redshift z. Following Natarajan, Croton & Bertone (20081, 



from the general formula the energy output integrated on 
the volume of the halo is given by 



< av > c 
2m^ 



(32) 



where poM is the complete DM density profile of the halo, 
accounting for the smooth main halo, the substructure and 
the adiabatic compression in both (see Section 3.2 - Section 
3.4). This luminosity can be decomposed into the main halo 
luminosity and the substructure luminosity such that: 



L 

< av > c 



2my 



/ pDM{r,Mmain,z)dV, 

Jv 



(34) 



and 



dN 
dm 



L^{m, z)dm. 



Note that in the above equations, L^{M, 



(35) 
is the 

DM annihilation luminosity of a smooth, adiabatically com- 
pressed halo or subahlo, the only difference among them is 
that for the halo or the subhalo, the mass of the SMBH de- 
pends differently on the mass of the halo. It is important 
to mention that in order to account for the substructure lu- 
minosity, an integration of the substructure mass function 
multiplied by the luminosity of an individual subhalo in all 
the substructure mass range must be performed. This inte- 
gration should be done between irifree, the free streaming 
mass of the WIMP (the minimum mass that a virialized DM 
halo can have), and Mmain, the mass of the main halo (no 
substructure can be as massive as the main halo). The free 
streaming mass will depend naturally on the WIMP mass, 
as mfree OC . Based on the analysis of the effects of the 
free streaming length on the matter power spectrum by |An-| 
gulo & White (20101, we consider a free streaming mass of 



lO-'^M, 



for a WIMP mass of = ITeV, and a 



free streaming mass of = 10 Mq for a WIMP mass of 
= lOGeV. 

We explicitly compute the luminosity of DM haloes, 
in the ACDM cosmology, considering two different WIMP 
masses (lOGeV and ITeV) and the presence or absence of 
substructure and adiabatic compression. We show the re- 
sults obtained in Figure 3. 

It can be seen that, for the same halo mass, the luminos- 
ity of the halo increases with redshift. This is a result of the 
structure formation process, because a DM halo that formed 
earlier will be more concentrated and have a higher char- 
acteristic density. Thus, its annihilation rate, and in con- 
sequence its luminosity, will be increased. Also, it can be 
seen that the greatest effect on the luminosity comes from 
the WIMP mass. In accordance with what was found for 
the smooth case, is almost inversely proportional to m^. 
There are variations to this proportionality, in the cases with 
substructure and adiabatic compression, because the maxi- 
mum attainable mass of the compressed spike, as well as the 
minimum mass (free streaming mass) of the substructure are 
also dependent on m^. These deviations are, however, min- 
imal. It can be also seen that the second greatest effect on 
the luminosity comes from the presence of adiabatic com- 
pression, and that the presence of substructure is the less 
important boosting mechanism. 

It is also important to mention that according to 



Natarajan, Croton & Bertone (20081, the energy output re- 



(33) 



quired to offset the cooling flows in galaxy clusters (in order 
to quench star formation at low redshifts in galaxy cluster 
systems, and recover the properties of the Red Sequence), 
that in current semi-analytic models of galaxy formation is 
provided by AGN feedback, is on the order of 10^^ - lO^V/s 



12 Ignacio J. Araya and Nelson D. Padilla 



10* 10'° 
M [Mj 



M [Hal 



10* 10* 10'° 10'° 

M [Mj 




M [MJ 



Figure 3. The internal DM annihilation luminosity of individual 
haloes, as a function of halo mass (in the range from 10""' to 10^'' 
Mq ) . The different curves in the four plots correspond to different 
redshifts; in particular, the dotted curves correspond to z = 0, the 
short-dashed curves correspond to z=3.6, the dot-dashed curves 
correspond to z = 10, the triple-dot-dashed curves to z = 50.1 and 
the long-dashed curves to z = 100. The different plots correspond 
to different clustering scenarios and WIMP masses; in particular, 
the top-left plot corresponds to the case with substructure and 
adiabatic compression, with a WIMP mass of ITeV; the top-right 
plot corresponds to the case with substructure and ITeV WIMP 
mass, but with no adiabatic compression; the bottom-left plot 
corresponds to the case with a WIMP mass of ITeV and adi- 
abatic compression, but no substructure; and the bottom-right 
plot corresponds to the case with adiabatic compression and sub- 
structure, but a WIMP mass of lOGeV. In all the plots we also 
include the v = \ mass values for the 3 lowest redshifts in the 
corresponding line-styles, and the 10''^J/s luminosity threshold 
corresponding to AGN feedback in galaxy cluster-sized systems. 



and thus, in general, the luminosities that we obtained are 
lower than the required value (by some 2 — 4 orders of mag- 
nitude, except in the case of = lOGel/, but this needs 
further study). It should be remembered that, although we 
plot halo luminosities up to haloes of mass of 10^^, the most 
massive gravitationally bound haloes today (at z = 0), have 
masses of the order of lO^M©, so more massive haloes 
are virtually non-existent; and further more, the masses of 
the haloes that have a particular number density at a certain 
redshift (given by the characteristic dimensionless mass (yY), 
falls rapidly with redshift. For clarity, we also include in Fig- 
ure 3, the lines corresponding %o v = \ sX different redshifts 
and the 10'^^ J/ s luminosity threshold of AGN feedback in 
galaxy cluster-sized systems. 



4.3 The global, clumped, DM annihilation 
luminosity per baryon 

Having computed the DM annihilation luminosity of a 
particular halo of mass M at redshift z, we now return to 
the problem of computing the global energy injection rate 
to the IGM per baryon, but considering the dumpiness of 
the DM on all clustering scales. To do this, we simply con- 
sider that, knowing the luminosity of a particular halo (that 



already includes substructure and adiabatic compression), 
and knowing the halo mass function (calculated in Section 
3.1), the energy output rate density (or luminosity density) 
due to DM annihilations is given simply by: 



clumped 



(Z) = (1 + ^)^ 



LAM, 



dn 



where m/ree is the free-streaming mass, 



{M,z)dM, (36) 



is the maxi- 



minosity (main halo plus substructure) of a DM halo of mass 



M at redshift z. 



is the previously calculated halo mass 



function (in comoving coordinates), and the factor {1 + z) is 
to convert the comoving density to a proper density. We con- 
sider nimax ~ W^"^ Mq since the number density of collapsed 
DM haloes with mass above this limiting mass is negligible. 

Finally, the energy injection rate to the IGM per baryon 
can be directly computed as: 



e: 



clumped 



clumped 



(^) 



nb(z) 



(37) 



where nb{z) is the proper number density of baryons at red- 
shift z (equal to (1 -I- z^^rihfi in terms of the present day 
baryon number density computed in Section 4.1). 

It is important to mention that the energy injection rate 
to the IGM per baryon computed in this way, considers only 
the contribution of the DM haloes and not of the smooth 
component. In the case of a zero free-streaming mass, all 
the DM would be in haloes since there would be no min- 
imum mass scale for the formation of collapsed structures. 
However, due to the non-zero free streaming mass, all haloes 
with masses less than that value, that should exist according 
to the halo mass function, are simply not allowed to form. 
As the clumped energy injection should be always greater 
than the smooth energy injection (computed in Section 4.1), 
and as the contribution due to DM haloes will be declining 
rapidly with increasing redshift after the truncation of low 
mass haloes due to the free-streaming limit becoming im- 
portant, a good approximation is to use as the clumped DM 
injection rate the greater value among the one calculated in 
equation (47), and the one calculated in Section 4.1. This is 
the approach that was used for the case of a IQGeV WIMP. 

The results obtained for the clumped energy injection 
rate per baryon, considering different clustering scenarios 
and WIMP masses, are shown in Figure 4. We only con- 
sider the maximal absorption fraction [fabs ~ 1), because for 
other values, the resulting curves should simply be rescaled 
by fabs- 

It can be seen that the case considering a WIMP with 
mass of lOGeV is the one that gives the higher injection rate, 
in agreement with what was discussed above, for redshifts 
lower than a certain value of z ~ 50, and gives the lower 
injection rate for higher redshifts up to z ~ 100. This can 
be understood as a consequence of the higher free-streaming 
mass (and thus, higher minimum mass that haloes can have), 
and the fact that the mass of the haloes that have a charac- 
teristic v = 1 (the haloes that have just recently collapsed) 
decreases with increasing redshift, such that with a higher 
free-streaming mass, less haloes are allowed to form at higher 
redshifts. Also, the presence of adiabatic compression is im- 
portant for boosting the energy injection rate, but only at 
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Figure 4. The clumped energy injection rate per baryon {E^{z)) 
as a function of redsliift, for different scenarios, in thick lines. The 
dotted curve corresponds to the case with substructure and adi- 
abatic compression for a WIMP mass of ITeV, the long-dashed 
curve corresponds to the case with ITeV WIMP mass and sub- 
structure but no adiabatic compression, the dot-dashed curve cor- 
responds to the case of ITeV WIMP with adiabatic compression 
but no substructure, and the short-dashed curve corresponds to a 
WIMP mass of lOGeV with substructure and adiabatic compres- 
sion. Also included as thin lines are the previous curves of Figure 
2, corresponding to the smooth cases, for a ITeV WIMP (dotted) 
and a lOGeV WIMP (short-dashed). 



Figure 5. The dumpiness factor C(z) as a function of redshift. 
The thick curves correspond to our results. The line-styles of the 
thick curves are given as in Figure 4. For comparison, we also 
include the curves obtained by Cumberbatch, Lattanzi & Silk 
(2010) for the case of an NFW profile, taken from the upper panel 
of their Figure 3, as thin curves. The line-styles of the thin curves 
correspond to different values of {Mmin/ModotMcut/Modot)- 
The solid line corresponds to (10^,10'^), the dotted line cor- 
responds to (10"'', 10®), the long-dashed line corresponds to 
(10-'',10-*), the short-dashed line corresponds to (10"^^, 10'^) 
and the dot-dashed line corresponds to (10~^^,10~'^). We refer 
the reader to the paper of the authors for further details. 



redshifts lower than 10. The presence of substructure is the 
least important of the boost mechanisms. 

In the next section, we obtain the dumpiness factor, 
C(z), implied by our energy injection calculations, and com- 
pare our results with those obtained by [Cumberbatch, Lat-| 
tanzi & SilklpoTol). 



4.4 The dumpiness factor 

The dumpiness factor is a useful tool when accounting 
for DM clustering in global DM annihilation energy injec- 
tion computations, like the ones required for calculating the 
heating and ionization of the IGM due to DM (see Section 
6.1 and Section 6.2). It simply relates the clumped energy 
injection to the smooth one by: 



C{z) = ^ ^. 

^smooth (2) 



(38) 



We show the values obtained for the dumpiness factor 
C(z) in Figure 5. The meaning of the curves are explained 
in the caption. 

In general, it can be seen that for all the cases consid- 
ered, the dumpiness is a decreasing function of redshift. This 
is reasonable, because as the redshift increases, the universe 



will tend to be more homogeneous and uniform, because 
the gravitational collapse and formation of structures would 
have had less and less time to occur. Also, for all the cases 
considered, the universe will be practically completely ho- 
mogeneous by redshift (2 ~ 150), except for the case with 
lOGeV DM WIMP, in which the universe becomes homoge- 
neous at redshift (z ~ 70). Note that in this case, although 
■^clumped ^^-^ is higher than for the other cases by almost an 



order of magnitude at redshift zero, i?J 



^(z) is higher 



as well, and so the dumpiness C(z) is not. It can also be 
seen that in agreement with what was discussed before, the 
adiabatic compression mechanism is important for low red- 
shifts, and the presence of substructure is the less important 
mechanism for increasing the dumpiness. 

We show also the dumpiness factors C{z) obtained by 
Cumberbatch, Lattanzi & Silk (2010), in Figure 5, using 



different line-styles. The two sets of curves can be readily 
compared as the definitions of the dumpiness factors are 
the same. There are, however, differences in the methods 
for computing them. In particular, [Cumberbatch, Lattanzi| 
[fc Silk[ ( |2010[ ) compute the halo mass function using the 
P-S theory, starting from an analytic fit to a{M,z), and 
then consider that substructure and sub-substructure follow 
a pure power-law mass function with exponent —2. We re- 
fer the reader to the authors for further detail about their 
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computation procedure. In general, it can be seen that we 
predict less dumpiness than the authors at low redshift, but 
higher dumpiness than the authors at high redshift. This 
can be understood because the authors also include sub- 
substructure (sub-subhalocs within subhalocs), that may be 
more important near redshift zero, but we include adiabatic 
compression (and they do not), which may be important at 
higher redshifts than substructure, but low enough so that 
formation of haloes is not disfavoured by the free-streaming 
mass. 

Having computed the energy output due to DM annihi- 
lations of local haloes, and having obtained the energy injec- 
tion rate into the IGM per baryon in the perfectly smooth 
case and in the clumped case, considering a standard ther- 
mal relic WIMP cross section, we consider in the next section 
the Sommerfeld enhancement, a physical process that may 
be responsible for enhancing the annihilation cross section 
by as much as 3 orders of magnitude. 

We do not include Sommerfeld enhancement in our 
calculations, because it is highly dependent on unknown 
physics (the mass of the new force carrier and the result- 
ing potential, as well as the physical scenarios in which this 
new force could arise), however, its effect is to directly mod- 
ify < av > hy & constant factor. As we already mentioned 
in the Introduction, the presence of Sommerfeld enhance- 
ment is controversial and not supported by all observational 
evidence. 

Having computed the DM annihilation luminosity gen- 
erated within a DM halo for different scenarios, in the next 
section we study the luminosity received by a DM halo due 
to the rest of the universe. 



5 EXTERNAL DM ANNIHILATION ENERGY 
INJECTION ON A HALO 

On previous sections, we iiavc considered particular 
WIMP models, calculated the DM clustering on all scales 
explicitly, and computed the DM annihilation luminosity 
generated within a particular DM halo, and the energy in- 
jection per baryon on the IGM. In this section, wo calculate 
the DM annihilation luminosity received by a halo due to the 
energy generated in the rest of the universe. For this, we will 
use the formalism developed so far, and calculate the energy 
output of a distant halo that is received by a particular halo. 
Also, we will use the bias factor and the two-point correla- 
tion function, that allow to calculate the spatial clustering 
statistics of haloes around a particular halo. Then, combin- 
ing these, we will calculate the desired energy received by 
the halo from all DM structures external to it. 

In the next subsection, we consider and compute the 
luminosity that is received by a distant halo. 



5.1 The luminosity received by a distant halo 

The flux received by an astronomical source, located at 
a nearby distance d, follows the usual relation: 



f- 



(39) 



should be replaced by the luminosity distance dL{r), that is 
a function of the comoving distance r to the source (equal 
to the proper distance at redshift zero). Also, if one wants 
to compute the received luminosity (energy per unit time) 
from the source, one should multiply the flux / by the cross 
section to photons of the surface that receives the radiation 
(e.g., ttR^ for a sphere). Finally, if the travelling photon has 
a non-zero probability of being absorbed in the way, and 
thus may not arrive at the receiving surface, this should be 
accounted for by the addition of a certain transfer function. 
Thus, the desired luminosity that is received by a distant 
halo can be written as: 



Lx,i{M,Mi,z,r) 



(M) 

47rdi(^,r) 



(40) 



where L-^,i{M,Mi,z,r) corresponds to the luminosity that 
is received by a halo of mass M at redshift z, due to a halo 
i of mass Mi located at a comoving distance r from the first 
halo. The function frir) corresponds to the transmittance 
function, (Tbary{M) corresponds to the baryonic cross-section 
of the receiving halo, Zi{z,r) corresponds to the redshift at 
which a photon is emitted by halo i such that it arrives to 
the receiving halo at redshift z, and dL{z,r) is the lumi- 
nosity distance between two halos separated by a comoving 
distance r at redshift z. 

For obtaining the redshift of emission {zi{z, r)), we con- 
sider that the comoving distance travelled by a photon emit- 
ted at redshift Z2 and received at zi is given by: 



/ 1+2=1 cda 
1 ^?H(a 



(41) 



This gives Zi — Zi{z, r) as an implicit function that can 
be evaluated numerically. 

For computing the luminosity distance between two ob- 
jects separated by a comoving distance of r at redshift z (not 
to be confused with the luminosity distance of an object at 
redshift z with respect to an observer at redshift zero), we 
use 



dz,(r, z) 



1 + z 



Sfe(r) 



{1 + Zi{z,r)) 
1 + z ' 



(42) 



where the function Sk{r) corresponds to the cosmologi- 
cal distance measure (given by the angular part of the 
Robertson- Walker metric), and is different for flat and 
curved cosmologies. For a flat case (Einstein-De Sitter and 
ACDM cosmologies), Sk{r) = r. 

For the halo cross section, we consider a baryonic cross 
section given by: 



(Tbary{M) = {rbary,Mw) T- 



JV^2/3 



M 



2/3 



(43) 



In the case of cosmological distances, the distance d 



where rbary,Mw = O.OlMpc and MMw,haio = 10 Mq are 
fiducial values for the baryonic radius and halo mass of a 
MW-like DM halo. We therefore consider that the cross sec- 
tion of photon absorption of a halo scales as its mass to the 
2/3 power, or in other words, as its transversal area. We also 
consider that the baryons are responsible for the photon ab- 
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sorption, that the cross section is redshift-independent, and 
that it does not depend on halo orientation. This cross sec- 
tion is in proper Mp(? . 

The transmittance function was obtained by analogy to 
radiative transfer theory. If the cross section of a DM halo 
depends only on its mass, we can define an optical depth for 
DM-annihilation-generated photons travelling through the 



universe as: 



r{z) 



a{M)^{l + zfdM, (44) 



where the (1 + z)^ factor is to convert the cross section to 
comoving coordinates. 

Then, the transmittance function is simply given in 
terms of the optical depth by: 



with masses above M, and ^DMir,z) the unbiased corre- 
lation function (labelled simply as ^{r,z), and considered 
previously). This bias factor is independent of r in a first 
order approx. 

We adopt the bias factor given by |Sheth, Mo fc Tormen| 



(20011 that takes into account ellipsoidal collapse and fits 



well numerical simulation measurements. 



5.4 The external luminosity received by a halo 

We now have all the pieces required to compute the 
external luminosity received by a DM halo due to all DM 
haloes external to it. First, we consider the radial halo 
mass function, that will be given by the mass function com- 
puted with P-S theory, but considering the clustering of DM 
around the particular halo of mass M: 



fTir,z) = e 



(45) 



where eabs is the fraction of photons that are indeed ab- 
sorbed when they pass through the cross-section of a halo. 

Now, all the functions used in the computation of L^^i 
(the luminosity received due to halo i) have been specified. 
In the next section, we explain the two-point correlation 
function that measures the spatial clustering of haloes. 



5.2 The two-point correlation function 

The correlation functions we will consider take into ac- 
count redshift evolution 

Although the correlation function is completely spec- 
ified by the power spectrum, the P{k) that we calculated 
is valid in the linear theory approximation, and thus, it is 
not a good estimator of the DM halo spatial clustering on 
small scales. Therefore, we will use a standard analytic fit to 
the zero redshift correlation function, as a pure power law, 
with its exponent and normalization motivated by observa- 
tion (for example the Sloan Digital Sky Survey, as analysed 
by Zehavi et al. (2011 1) and simulations: 



Co(r) = 



(46) 



where r is in comoving Mpc. 

In the next section we explain the bias factor, that rep- 
resents the excess in clustering of haloes above different mass 
scales. 



5.3 The bias factor 

In the previous section, we defined the correlation func- 
tion ^(r, z) as the excess clustering of the matter density 
field on characteristic separations r at redshift z. DM haloes 
with different masses (above a mass scale M), provide biased 
measurements of clustering. The bias factor of mass scale M 
can be defined by: 



b^{M,z) 



£,MM (r, z) 



(47) 



£,DM{r-, z) ' 

with £,MM{r,z) the correlation function of the DM haloes 



dn 



{M,M„z,r) = {l + b{u{A-U,z))b(u{M,z)) x 

dn 



Now we define the radial emissivity, as the DM anni- 
hilation energy generation rate per comoving volume at a 
distance r from the halo of mass M: 



e^{M,z,r) 



/1 71 

J m 



dn 
dMi 



{M,Mi,z,r) X 



xLxAM,Mi,z,r)dMi. (49) 
Finally, the required external luminosity will be given 



by: 



(M, z) = eabsEe 



4TvSl{r)e^{M,z,r)dr. (50) 



The volume integration is performed between the co- 
moving virial radius of the halo up to rmax ■ Formally, r^ax 
should be oo in the Open and Flat cosmologies, and ttRq in 
closed cosmologies. However, in practice, we use a radius of 
rmax = IGpc, because at such large distances, the luminos- 
ity of individual haloes arrives so diluted that the contribu- 
tion to the integral is negligible. We also consider the ab- 
sorption and escape fractions eats and eesc- This absorption 
and escape fractions are different to the previously discussed 
fabs (z) , because they are the fraction of the DM annihilation 
energy that is trapped (absorbed) within haloes (eats), and 
the fraction of the DM annihilation energy that manages 
to escape from haloes as photons and is free to propagate 
through the universe. We assume that eats + ^esc = 1, so that 
no energy is wasted, and consider eabs ~ 0.9 and eesc = 0.1. 
We choose these values because an escape fraction of 0.1 
means that 10% of the photons manage to reach the IGM, 
and therefore, this choice is somewhat similar to choosing 
fabs{z) = 0.1 as motivated by Cumberbatch, Lattanzi & 



Silk (20101 and Natarajan fc Schwarz (20091, although both 



fractions are not directly comparable. 

In Figure 6, we show our results for the external DM 
annihilation luminosity as a function of halo mass at dif- 
ferent redshifts, considering only the ACDM cosmology, in 
the case of a WIMP with mass ITeV and considering both 
substructure and adiabatic compression. 
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to their extra energy injection (for example, see |Pierpaoli 


(20041; Mapelli, Ferrara & Pierpaoli 


(20061; Ripamonti, 


Mapelli & Ferrara (20071; Valdes et al. 


(20071; Natarajan & 


Schwarz| (|2009|); and |Cumberbatch, Lattanzi & Silk| (|2010|) 



Figure 6. The external DM annihilation luminosity received by 
individual haloes, as a function of halo mass (in the range from 
10~* to lO^'^ Af©)- The colours of the curves correspond to dif- 
ferent redshifts; in particular, the dotted curve correspond to red- 
shift zero, the short-dashed curve correspond to redshift 3.6, the 
dot-dashed curve correspond to redshift 10, the triple-dot-dashed 
curve to redshift 50.2 and the long-dashed curve to redshift 100. 
The curves were computed assuming a WIMP mass of ITeV and 
considering both substructure and adiabatic compression. We also 
include the v = \ mass values for the 3 lowest redshifts in the 
corresponding line-styles, and the 10'^^ J/s luminosity threshold 
corresponding to AGN feedback in galaxy cluster-sized systems, 
as in Figure 3. 



In general, comparing these results with the ones ob- 
tained for the DM luminosity generated within a particu- 
lar halo (see Figure 3), one can see that the external lumi- 
nosity relation is less steep than the internal one, and so, 
therefore, the external luminosity will be more important 
than the internal luminosity for haloes below a certain mass 
(~ 10'^ — lO^M©), and less important for haloes above that 
mass. The implications of this have not yet been studied, 
but it will be interesting to do so as future work. 



6 GLOBAL IGM HEATING, IONIZATION 
AND 21 CM SIGNAL DUE TO DM 

Different authors have studied the effects of DM annihi- 
lations (or decays) on the global properties of the IGM, due 



among others). In this section, we study the effects on the 
IGM of the energy injection found on Section 4.3. We divide 
the section into different subsections, each corresponding to 
one of the different physical quantities of interest that are 
affected by this extra energy injection. We explain the for- 
malism and the definition of the quantities, as well as how to 
compute them, in the beginning of the corresponding sub- 
sections, and then we show the results we obtained and com- 
pare them with different authors. The first two sections (the 
ionization fraction and the IGM temperature), refer to gen- 



eral properties of the IGM, whereas the next sections have 
to do with the signal of this extra- heating on the 21 cm line 
of the hyperfine structure of neutral hydrogen (HI). 

We start, in the next subsection, with the computation 
of the IGM ionization fraction. 



6.1 The IGM ionization fraction 

The ionization fraction (a;e) is defined as the ratio of 



free electrons to hydrogen atoms a;^ = 



Following 



Valdes 



et al. (20071, the variation of the ionization fraction is given 



(51) 



by: 



dXe. _ 

~ H{z){l + z) 

where Is and Rs are the standard ionization and recombina- 
tion rates per baryon, and = Xi i^) ^ is the contribution 
to the ionization rate due to the energy injection of DM an- 
nihilation, is the energy injection rate per baryon com- 
puted in Section 4.3, Eo is the ionization energy of Hydrogen 
(13.6 eV), and Xii^) is the ionization efficiency and is given 
by Xii'^) = ^~3^° (i-e., if all matter is ionized, no energy is 
used for ionization, and if all matter is neutral, 1/3 of the 
injected energy goes to ionizations). The factor u{z){\+z) ' 
is simply a change of variables, because the differential of 
cosmic time dt and the redshift differential dz are related by 
dt ■ 



-dz 



- H{z)(l + z)- 

We define the DM annihilation contribution to the ion- 
ization fraction as 5xe{z) = Xe{z) — x"°~^^^ (z), where 
xr-°" is the ionization fraction in the standard scenario. 
It is important to mention that this standard ionization frac- 
tion is a consequence of the physics of recombination from 
the epoch of the last scattering surface, and is simply the 
relic free-electron abundance from the freeze-out of the re- 
combination reaction. Thus, the j:""^^*^ does not account 
for the reionization of the universe by Poulation HI stars, 
that according to WMAP7 ( Komatsu et al.^2 011), occurred 
at redshift 10.6 ± 1.2. With this definition of Sxe, we have: 



dSx^ 
dz 



-Jx(^) 
Hiz)il + z) 



(52) 



We solve this differential equation, using the x"° 
given in Figure 3 of |Mapelli , Ferrara & Pierpaoli (20061, 



by the Euler method, using as boundary condition that 
5xs{z — 1000) — (at redshift z > 1000, as we are near 
the recombination epoch, the ionization fraction is near 1 
and is almost entirely due to the standard cooling scenario) . 

We show the obtained ionization fraction (xeiz)) curves 
in Figure 7 and in Figure 8. The first figure shows the Xe{z) 
curves obtained for the different DM clustering scenarios 
and WIMP masses, but considering the clumped DM case, 
and with an absorption fraction of fabs = 1 (more details 
on the caption of the figure). The second figure shows the 
ionization fractions, assuming a WIMP mass of ITeV and 
substructure and adiabatic compression, but for different 
values of the absorption fraction {fabs = 0.1 or 0.01), or 
simply not considering DM clustering at all, and using the 
homogeneous DM energy injection computed previously. 

It can be seen from Figure 7 that the ionization fraction 
evolves very differently from the case with no DM energy in- 
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Figure 7. The ionization fraction of the IGM for different cases, 
assuming a maximal absorption fraction and the dumpiness of the 
DM. Our results are shown as thick lines. The different line-styles 
of the thick lines correspond to: (short-dashed) a WIMP mass of 
lOGeV with adiabatic compression and substructure, (dotted) a 
WIMP mass of ITeV with adiabatic compression and substruc- 
ture, (long-dashed) a WIMP mass of ITeV with substructure but 
no adiabatic compression, (dot-dashed) a WIMP mass of ITeV 
with adiabatic compression but no substructure, and (solid) the 
standard cooling scenario, with no DM energy injection. For com- 
parison, we also include the results obtained by Mapelli, Ferrara 
& Pierpaoli (2006) for the case of the lOOGeV SUSY neutralino, 
taken from the bottom panel of Figure 3 of the authors and 
we show them as thin lines. The thin solid line corresponds to 
< av >= 2 X 10~^^cm^s~^ and the thin dotted line corresponds 
to < av >= 10~^*cm^s~^. Also included is the line correspond- 
ing to ^ = 20, that indicates when the effects of Population III 
stars may become important. 



Figure 8. The ionization fraction of the IGM for different cases, 
assuming a WIMP mass of ITeV and substructure and clustering, 
or no DM dumpiness at all. Our results are shown as thick lines. 
The different line-styles of the thick lines correspond to: (dot- 
ted) the standard clustering scenario considering both substruc- 
ture and adiabatic compression with a ITeV WIMP and fats = 
f, (short-dashed) the same clustering scenario but considering 
fabs = O-li (long-dashed) the same clustering scenario but with 
fabs = 0.01, (short-dot-dashed) the same WIMP mass but consid- 
ering no DM dumpiness, (long-dot-dashed) a lOGeV WIMP with 
no DM dumpiness, and (solid) the standard cooling scenario, with 
no DM energy injection. For comparison, we also include the re- 
sults obtained by Mapelli, Ferrara & Pierpaoli (2006) for the case 
of the fOOGeV SUSY neutralino, taken from the bottom panel 
of Figure 3 of the authors and we show them as thin lines. The 
thin solid line corresponds to < av >= 2 X 10~^^cm^s~^ and 
the thin dotted line corresponds to < av >= 10~^*cm^s~^. Also 
included is the line corresponding to 2: = 20, that indicates when 
the effects of Population III stars may become important. 



jection, as long as the DM dumpiness is considered. Also, 
for the case of a WIMP with a mass of lOGeV, the ioniza- 
tion fraction can be increased up to Xe ~ 0.1 at redshift zero, 
whereas for the other cases that consider a ITeV WIMP, this 
partial re-ionization can only bring Xe to ~ 0.01 at redshift 
zero, hence leading to a difference of one order of magnitude 
of the final ionization fraction (note, however, that as the 
computation of the ionization fraction does not account for 
the Population III stars, that are thought to have re-ionized 
the universe at redshift z ~ 10, this results are only reli- 
able down to that redshift at most; however we use a more 
conservative limit of z = 20 in the plots, to account for 
the fact that the formation of the first stars started earlier 
than the actual redshift of reionization). It can also be seen 
that both substructure and adiabatic compression are con- 
siderably less important in determining the evolution of the 
ionization fraction, and that adiabatic compression is more 
relevant near redshift zero, whereas substructure is more im- 
portant between redshifts ~ 7— ~ 100. 



From Figure 8, it can be seen that even without consid- 
ering DM dumpiness, the case with a lOGeV WIMP consid- 
erably modifies the ionization fraction of the universe. It has 
the effect of raising the relic ionization at redshift (with- 
out accounting for stellar populations) by a factor of ~ 20. 
Also, the evolution of the ionization fraction is particularly 
sensitive to the absorption fraction fabs, such that reducing 
fabs from 1 to 0.1 results in a completely different evolution 
(and in a ionization fraction at z=0 an order of magnitude 
lower). The cases with fabs =0.1 and no dumpiness give 
evolutions of the ionization fraction that are very similar to 
the case were no DM energy injection is present. 

6.2 The IGM kinetic temperature 

We now consider the kinetic temperature of the IGM. 
Following [Valdes et al.| ( [2007[ ), the change in this tempera- 
ture is given by: 
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az 



H(z){l + Y + x„) 



3hH{z)(l + Y + x. 



(53) 



where Ij = ^^'^^(3m"^e)""°^ analogue to a thermal con- 

ductivity between the CMB and the IGM, Xh is the heating 
efficiency and is given by Xh ~ "'""'"3^'' and Y is the helium 
fraction by mass (taken as Y = 0.26). 



It is expected (see for example Cumberbatch, Lattanzi 



& Silk (2010l) that the IGM kinetic temperature tracks the 



CMB temperature at high redshifts, so we used the CMB 
temperature at redshift z = 1000 as the boundary condition 
for the IGM kinetic temperature, and solved the differential 
equation by the Euler method. The interpretation of this 
cooling equation is straightforward, as the first term on the 
right hand side corresponds to the adiabatic cooling of the 
IGM in absence of heat exchange or energy injection, and 
proceeds as T oc (1 + 2)^. The second term corresponds 
to the heat exchange between the CMB and the IGM, and 
the third term allows for the extra heating due to the DM 
annihilation energy injection. 

We show the obtained IGM kinetic temperature curves 
in Figure 9 and in Figure 10. The first figure shows the tem- 
peratures obtained for the different DM clustering scenarios 
and WIMP masses, but considering the clumped DM case, 
and with an absorption fraction of fats ~ 1 (more details on 
the caption of the figure) . The second figure shows the tem- 
peratures, assuming a WIMP mass of ITeV and substruc- 
ture and adiabatic compression, but for different values of 
the absorption fraction {fabs = 0.1 or 0.01), or simply not 
considering DM clustering at all, and using the homogeneous 
DM energy injection computed previously. 

In can be seen on Figure 9 that for all the cases that 
consider the dumpiness of the DM (assuming fats = 1), the 
cooling is very different from the standard adiabatic sce- 
nario. Significant re-heating occurs in the case of a lOGeV 
WIMP, where the IGM may be heated even to tempera- 
tures of lOOOif by redshift zero (although, remember that 
this cooling scenario does not account for energy injection 
due to stellar components and so it is reliable only down to 
redshift z ~ 20). Also, adiabatic compression is more impor- 
tant for the cooling history than substructure, and the extra 
heating due to this compression can amount to an order of 
magnitude of difference by redshift zero. 

Also, it can be seen on Figure 10 that the absorption 
fraction is critical for determining the particular cooling his- 
tory, and that the low-redshift heating due to DM annihi- 
lations can vary enormously with different choices for fabs- 
It can be seen that, in general, in the cases considering adi- 
abatic compression, the IGM kinetic temperature reaches a 
quasi-stationary state, and the final temperature attained is 
almost proportional to the absorption fraction. 

We also include the results for the ionization fraction 
and IGM kinetic temperature obtained by |Mapelli, Ferrara| 
& Pierpaoli ( 2006 1 , in order to compare them with our re- 



sults. 

It can be seen that for the cases with no DM dumpiness, 
we reproduce the results obtained by the authors. Also, the 
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Figure 9. The kinetic temperature of the IGM for different cases, 
assuming a maximal absorption fraction and the dumpiness of 
the DM. For comparison, we also include the results of Mapelli, 
Ferrara & Pierpaoli (2006), taken from the top panel of Figure 3 
of the authors. The different line-styles are assigned in the same 
way as in Figure 7. The z = 20 line is also included. 



boost they considered in the cross section (for example, ac- 
counting for Sommerfeld enhancement), has almost exactly 
the same effect than lowering the mass of the WIMP par- 
ticle. Furthermore, we find significant effects not accounted 
for by the authors when considering the dumpiness of the 
DM on all its levels. 



6.3 The neutral hydrogen spin temperature 

We now study the 21 cm signal from neutral hydrogen, 
and how it is modified by DM annihilation. We foll ow |Cum-| 
berbatch, Lattanzi & Silkl (|2010|) andlValdes et al. I (120071) in 



the explanation of the physics and the processes involved, 
and refer the reader to the above papers for further details. 

The 21 cm line corresponds to the transition energy be- 
tween the singlet and triplet hyperfine levels (n = 0, j = 
and n — 0, j — 1) of the hydrogen atom. The spin temper- 
ature of the neutral hydrogen is defined as the temperature 
that relates the populations of singlet and triplet states: 



ni 
no 



3e" 



(54) 



where Ts is the spin temperature, T* is the equivalent tem- 
perature of the 21 cm transition (corresponding to T, — 
where A = 21cm), ni is the triplet state occupation 
number, no is the singlet state occupation number and the 
factor of 3 is the degeneracy of the triplet state (i.e., 3). 

In general, the spin temperature is a weighted average 
of the IGM kinetic temperature and the CMB temperature. 
In absence of additional heating mechanisms, both the IGM 
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Figure 10. The kinetic temperature of the IGM for different 
cases, assuming a WIMP mass of ITeV and substructure and 
clustering, or no DM dumpiness at all. For comparison, we also 
include the results of Mapelli, Ferrara & Pierpaoli (2006), taken 
from the top panel of Figure 3 of the authors. The different line- 
styles are assigned in the same way as in Figure 8. The 2 = 20 
line is also included. 



kinetic temperature and the spin temperature are coupled 
to the CMB at high redshifts [z > 300), whereas at lower 
redshifts, the IGM cools faster than the CMB. Thus, the 
spin temperature will be lower than the CMB and the IGM 
will be visible in 21 cm absorption. 

There are two main mechanisms that couple the spin 
temperature and the kinetic temperature: Spin exchange col- 
lisions and Ly-alpha pumping (also called the Wouthuysen- 
Field (W-F) effect ( Hirata| [2006)). Spin exchange collisions 
correspond to the collisions of neutral hydrogen (responsible 
for the 21 cm line) with other neutral hydrogen atoms, with 
electrons or with protons. These collisions will tend to ther- 
malize the populations of both levels (singlet and triplet) 
to the kinetic temperature of the IGM, and thus will cou- 
ple Ts with Tfe. Lyman-alpha pumping, on the other hand, 
corresponds to the mixing (and thus thermalisation) of the 
singlet and triplet levels by Lyman-alpha transitions (from 
energy levels 1 to 2 and vice- versa) to the 2p level. Therefore, 
the mixing of hyperfine states through the W-F effect will 
depend on the intensity of the Lyman-alpha photon field, 
which in turn will depend on the kinetic temperature of the 
IGM and on the direct energy injection due to DM WIMPS. 
Thus, the W-F will also couple Ts with T^. 

At low redshift {z < 50), the IGM density is sufficiently 
rarefied so as to make the spin exchange collisions too im- 
probable, and thus Ts will tend to follow the temperature of 
the CMB again. 

As the mechanisms coupling Ts and TkicM depend on 
the ionization fraction Xe, and the DM energy injection rate 
per baryon , it is expected that the presence of an annihi- 



lating WIMP will alter the evolution of Ts with redshift, and 
thus, the detection of this signal may help to discriminate 
between different WIMP scenarios. 

The spin temperature at redshift z can be computed as: 



Ts = 



TcMB + VaTk + VcTk 



(55) 



where ya is the W-F coupling, and yc is the coUisional cou- 
pling to the IGM temperature. The coupling coefficients can 
be computed as: 



AioTk ' 



and 



Vc = 



CioT, 
AinTk 



(56) 



(57) 



where Aio — 2.85 x 10~^^s~^ is the spontaneous emission 
coefficient for the transition, and Pio and Cio are the radia- 
tive and collisional transition coefficients respectively (they 
are rates, with units of [s^^]). The collisional transition co- 
efficient is given by: 



ClO = kloTlHI + ne"ie 



(58) 



where uhi and rig are the neutral hydrogen number den- 
sity and the free electron number density respectively, and 
fciQ and 7e are the specific (one particle) transition coeffi- 
cients for collisions with neutral hydrogen atoms and elec- 
trons respectively. We have neglected the contribution due 
to collisions with protons, following |Cumberbatch, Lattanzi] 
& Silk (20101. The number densities rie and nni are easily 



obtained in terms of the above-defined ionization fraction 
Xe, the hydrogen and helium abundances by number [fn 
and fne) and the baryonic density of the universe nj, also 
defined previously. We use the following formulas for the 
coefficients fcio and 7^: 



kw = 3.1 X 10""rfe ='"ea;p(-32_R'/rfc), [cm^s" 
and 

log{-^e/crrC' s 

= -Qmi + QMog{Tk/lK)exp{-[log{Tk/lK)\'^-'' 



(59) 



/1800), 



IK <Tk< W^K, 



(60) 



log{'ye/cm^s "*") 



-9.607 + 0.5Zo<?(rfe/lA'),rfe < IK, (61) 



7.(Tfe > W*K)^'y4lO^K). 



(62) 



The fit to kio is given by|Kuhlen, Madau fc Montgomery] 



(20061 and the fit to je is given by Liszt (20011 
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We now calculate Pio as: 



16 TT JaCTa 

27 hvc ' 



(63) 



where Jc, is the isotropic specific intensity of the radiation 
field in the Ly-alpha line, Oa is the Ly-alpha transition cross- 
section in the monochromatic approximation, and Va is the 
frequency of the Ly-alpha transition photon. According to 



Laor (19961, the a a monochromatic cross-section is given 
by: 



Tre ni 

0"a — /l2./se , , 

TTleC ni + 712 



(64) 



where /se is the stimulated emission correction and is given 
by: 



where r is the optical depth of the neutral IGM at 21(l+z) 
cm, and can be calculated as: 



-nHi(z), 



(70) 



where uo = 1420MHz is the frequency corresponding to the 
21 cm transition in the rest frame of the source. The differ- 
ential brightness temperature is important because it gives 
a difference in brightness between the CMB and the neu- 
tral hydrogen averaged in the beam of the instrument, and 
will be given by the above formula when the beam is large 
enough that the signal equals the average global one. Thus, 
this is the quantity that can be measured observationally 
in the low frequency range, and will be used for the next 
section. 



U = (1 



(65) 



the first and second energy levels population ratio ^ of the 
hydrogen atom is given by: 



6.5 The difference of differential brightness 
temperatures 

Finally, we compute the difference in differential bright- 
ness temperature, defined as: 



n2 
m 



52 

— ( 
51 



(66) 



the n — th level degeneration is given by Qn — 2n , and /12 is 
the oscillator strength of the transition (/12 = 0.416). Note 



that and the factor e in should be replaced by 



when 



using MKS units. 

Finally, the specific intensity of the Ly-alpha radiation 
field is given by: 



Ja{z) 



AtvH{z) 



<iff 

^p^2^ P 



XcE^iz) 



nuhva 



(67) 

where Xp — (1 — xhi) is the proton fraction with respect 
to the total hydrogen number density, ct^ip is the radiative 
recombination coefficient to the n = 2, I = 1 level, given 



in Table 1 ofjPengelly (19641, 7e_ff is a coUisional excitation 
rate of HI atoms involving electrons, and Xa = where 
Xe is the excitations efficiency given by Xe = ^—f^ (here, it 
can be seen that Xe + + Xft — li s-nd thus, all the en- 
ergy injected to the IGM by DM annihilation goes to either 
heating, ionization or excitation). For 7eir, we use: 

7eH = 2.2 X lO'^exp [-11.84/(T/10''7^)] , [cm^s"^], (68) 



given by [Cumberbatch, Lattanzi fc Silk| ( |2010[ ) . 

Now, we can explicitly compute the spin temperature 
Ts as a function of z, and use it in the following calculations. 



6.4 The 21 cm differential brightness temperature 

Having obtained the spin temperature Ts, we proceed 
to compute the 21 cm differential brightness temperature 
STt, defined as: 



STt 



Ts — Tch 
1 + z 



(69) 



ASTi{z) = \STbiz) ~ 5Tt,o{z)\, 



(71) 



where 5Tb,o{z) corresponds to the differential brightness 
temperature obtained without the extra energy injection by 
DM annihilations. According to |Valdes et al.| ( |2007| ), this 
signal could be measured with current and future radio tele- 
scopes, such as LOFAR, 21CMA, MWA and SKA, in prin- 
ciple, up to the « ImK level at redshifts where baryonic 
energy injection processes (such as Population III star for- 
mation) are not important {z >~ 20 — 30). However, the 
authors explain that the various foregrounds (i.e., Galactic 
free-free and synchrotron emission, unresolved extragalactic 
radio sources, free-free emission from ionizing sources, syn- 
chrotron emission from cluster radio haloes and relics) are 
much stronger than the cosmological signal, and are also 
very difficult to remove. The authors also explain some pos- 
sible methods for effectively removing the foregrounds and 
for clearly identifying the extra energy injection signal; how- 
ever, that is beyond the scope of this work and we refer the 
reader to their paper. 

In Figure 11 and in Figure 12, we show the results we 
obtained for the ASTb{z) curve, for the different cases we 
considered. In Figure 11 we assume the maximal absorption 
fraction fabs = 1, and consider the different DM clustering 
scenarios and WIMP masses. In Figure 12, we consider the 
standard scenario for the DM clustering (adiabatic compres- 
sion with substructure), and we vary the absorption fraction, 
the WIMP mass, or consider the perfectly smooth case (for 
more details about the plots, see their captions). 

It can be seen, from Figure 11, that the most impor- 
tant factor that determines the ASTt signal is the WIMP 
mass. For a WGeV WIMP, the signal is around ~ iOlmK] 
at redshift 100, very much above the ~ ImK sensitivity 
of current and future instruments quoted by [Valdes et al.| 
( |2007[ ). However, the authors also mention that current and 
future radio frequency arrays will probe up to a minimum 
frequency of 30MHz, thus corresponding to a 21cm emis- 
sion at 2 ~ 46. Therefore, the sensitivity window should be 
considered as ASTt > ImK, and also ~ 20 < z <~ 46. The 
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Figure 11. The difFerence of 21cm differential brightness tem- 
perature {A5Ti,{z)) for different cases, assuming a maximal ab- 
sorption fraction and the dumpiness of the DM. Our results are 
shown as thick lines. The different line-styles are assigned in the 
same way as in Figure 7. For comparison, we also include the 
results obtained by Cumberbatch, Lattanzi & Silk (2010) for the 
case of SUSY models with a NFW density profile, taken from the 
upper panel of Figure 10 of the authors, and we show them as 
thin lines. The line-styles for the thin curves correspond to dot- 
dashed for Model 4, short-dashed for Model 2, long-dashed for 
Model 3 and dotted for Modell. We refer the reader to Table I 
of the authors for a description of the different models. We also 
include the detectability limits here discussed, corresponding to 
the 1 mK threshold and to the z = 20 (Population III stars) and 
z = A6 (minimum frequency of 30 MHz) limits. 



Figure 12. The difference of 21cm differential brightness tem- 
perature {A5Ti,{z)) for different cases, assuming a WIMP mass 
of ITeV and substructure and clustering, or no DM dumpiness at 
all. Our results are shown as thick lines. The different line-styles 
are assigned in the same way as in Figure 8. For comparison, 
we also include the results obtained by Cumberbatch, Lattanzi & 
Silk (2010) for the case of SUSY models with a NFW density pro- 
file, taken from the upper panel of Figure 10 of the authors, and 
we show them as thin lines. The line-styles for the thin curves 
correspond to dot-dashed for Model 4, short-dashed for Model 
2, long-dashed for Model 3 and dotted for Modell. We refer the 
reader to Table I of the authors for a description of the different 
models. We also include the detectability limits here discussed, 
corresponding to the 1 mK threshold and to the z = 20 (Pop- 
ulation III stars) and 2 = 46 (minimum frequency of 30 MHz) 
limits. 



IQGeV case is well within this window, due to its broad, 
secondary peak around z ~ 30, caused by the DM cluster- 
ing. This same feature can be seen in the AJTb curve of the 
ITeV case with adiabatic compression and substructure, al- 
though the intensity of the signal is considerably lower, and 
the peak is now the primary peak. In the case of the ASTb 
signal, the presence of substructure is more important than 
the presence of adiabatic compression. At least in principle, 
this signal could be detected for all the cases with fats = 1 
that consider dumpiness. 

It can be seen in Figure 12, that even in the no- 
clumpiness case, a WIMP with mass of lOGeV would still 
be inside the sensitivity window mentioned before, due to 
the broad peak in the signal around z ~ 100 and the sec- 
ondary peak around z ~ 25. Also, for the cases with ITeV 
WIMP but non-maximal fats or no dumpiness, the signal 
would be much weaker (on the ~ 1 — 3mK level) and would 
just barely fall in the sensitivity window (at 40 < z < 46), 
because the diminishing of the absorption fraction or the 
dumpiness would tend to shift the signal to redshift ~ 100. 

It is important to keep in mind that, as was men- 
tioned above, the sensitivity window here considered is only 



nominal, because the foreground extraction and the de- 
contamination of the signal may be difficult. 

We also present in Figure 11 and Figure 12 the A5Ti,{z) 
curves obtained by |Cumberbatch, Lattanzi fc Silk, ( |2010| 
and compare them with the results that we obtain. 

It can be seen from the comparison between the curves 
given by the authors and ours, that the dynamical range 
of the values for the obtained ASTt{z) curves are similar, 
but that the particular shapes of the curves are different. 
This may be because the authors use an analytic approx- 
imation to a{M,Z), and assume a pure power-law shape 
for the substructure, without considering time evolution of 
substructures within particular haloes, and so their curves 
will be more soft, having less numerical noise from the in- 
tegration and equation solving procedures. For the case of 
maximal absorption fraction, and standard clustering, with 
a WIMP mass of ITeV, we find that the signal is peaked 
around redshift 50, in agreement with the results obtained 
by the authors, but for the cases with non-maximal absorp- 
tion and less or no clustering, our curves are more displaced 
towards z~100. It is interesting to note, however, that the 
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curves presented by the authors aheady consider the absorp- 
tion fraction explicitly. 



7 CONCLUSIONS 

In this paper, we develop formalisms for calculating the 
generated and received DM self annihilation luminosity of 
individual halos, as well as the energy injection rate per 
baryon in the universe as a whole, accounting for the clus- 
tering of the DM at all scales but without considering Som- 
merfeld enhancement. We also compute the effects on the 
IGM and on the 21cm signal of this extra energy injection 
due to the DM. 

It is important to mention that we take a numeric 
approach for performing the calculations of the barionic 
dumpiness, instead of the analytic approach used by pre- 
vious authors. Also, we incorporate adiabatic compression, 
which had not been done so far. The prescriptions used 
for the computations are the primordial power spectrum, 
the transfer function, the window function, the substruc- 
ture mass function and its corresponding time evolution, the 
computation of the NFW parameters as a function of halo 
mass and redshift and the computation of the adiabatically 
compressed density profile for each halo. 

We now give some general conclusions about the results 
obtained. 

As discussed in Section 4.2, the luminosity generated 
within individual halos, considering ITeV WIMPs, is not 
enough to self-regulate the star formation of the halo, and 
thus, it cannot contribute significantly to the standard feed- 
back mechanisms (AGN feedback for galaxy cluster-sized 
haloes of mass greater than lO^^M©). In the case of a 
lOGeV WIMP and considering clustering by substructure 
and adiabatic compression, for haloes of masses greater than 
IO'^^Mq, dm annihilation may indeed contribute signifi- 
cantly (even more than AGN feedback) to the total energy 
injection into the Intra-Cluster Medium (ICM). The latter 
case however, considers a WIMP mass that is not the most 
favoured one. Also, virialized DM haloes of masses greater 
than IO^^Mq are rare even at redshift zero. At higher red- 
shifts, high mass haloes are increasingly rare, and therefore, 
the importance of DM annihilation as a feedback mechanism 
is further diminished. 

DM clustering in general is determinant in order to as- 
sess the effects of DM annihilation on the IGM and for com- 
puting the energy injection per baryon. Evidence of this are 
the dumpiness factors obtained in Section 4.4, which can be 
m excess of 10*^ at redshift zero (the energy injection rate is 
thus boosted by this factor with respect to the case of a per- 
fectly smooth universe). Also, in favour of the importance of 
clustering are the results obtained for the ionization fraction 
Xe{z) and the IGM kinetic temperature TkiGM{z) discussed 
in Section 6.1 and Section 6.2, where it is clearly seen that, 
for the cases assuming a perfectly homogeneous universe, 
and for a WIMP mass of ITeV, the obtained evolutions 
of these quantities are virtually indistinguishable from the 
standard case with no energy injection due to DM. For ex- 
ample, not accounting for substructure, led |Mapelli, Ferrara| 
&i Pierpaoli (20061 to conclude that the effects of a heavy 



this may not be the case by accounting for the clustering of 
DM. 

We find that the presence of annihilating and clumped 
DM may result in significant deviations in the evolutions of 
the inferred temperature and ionization fraction. For exam- 
ple, in the scenario that takes into account the clustering 
of DM at all levels, as well as the adiabatic compression, a 
ITeV WIMP with thermal relic annihilation cross-section 
may, for example maintain the temperature of the IGM on 
the ~ WOK level, or contribute to the ionization fraction 
maintaining it at a value of ~ 10~^, until redshift zero (not 
considering POPIII stars). We conclude that although the 
DM annihilation can have a significant impact on the in- 
ferred properties of the IGM, it cannot be regarded as an 
alternative reionization scenario, because even in the case of 
a lOGeV WIMP (with lower mass than favoured by indirect 
detections and accelerator constraints), the value of Xe at 
z < 100 does not rise over 0.1. This is in agreement with 



the results obtained by Mapelli, Ferrara & Pierpaoli (20061, 



(~ lOOGeV^) neutralino on the IGM were negligible, and it 



was shown by Cumberbatch, Lattanzi & Silk (20101 that 



where for all the different models of decaying or annihilating 
DM considered by the authors, they obtained values of Xe 
at 2 < 100 of at most 0.3 (for the case of axino-type Light 
Dark Matter with a mass of lOMeV and without considering 
POPIII stars). We note however that among the models that 
they considered, only the neutralino DM model is compara- 
ble to our results for the no dumpiness case. Both results are 
in good agreement, however we do not use exactly the same 
mass as they do, and we do not consider cross-section en- 
hancement. Varying the mass of the WIMP and considering 
a constant enhancement in cross-section can have enormous 
effects on the properties of the IGM as the DM energy in- 
jection rate scales as oc ^^"^ (disregarding the variation 
in the presence of substructure when changing due to 
free-streaming effects) . 

About the detectability of the 21cm ASTt signal, de- 
fined as the absolute value of the difference of the differential 
brightness temperatures in 21cm with and without DM en- 
ergy injection (see Section 6.5), we can conclude considering 
the results that we obtain that, according to the sensitiv- 
ity limits mentioned by [Valdes et al.| ( |2007[ ) , and consider- 
ing a maximal absorption fraction fabs = 1, a thermal relic 
WIMP with mass of ITeV may have a signal detectable at 
the ~ 15mA' level, peaked around z ~ 50 (corresponding 
to a frequency of ~ 30MHz; see Section 6.5). For lower ab- 
sorption fractions, a ITeV WIMP may not be detectable by 
the 21cm signal. A lower mass WIMP may be detectable, 
although a mass on the order of ~ lOGeV^ or lower is not 
favoured by neither accelerator constraints (see Section 2.3), 
nor indirect detection searches (see Section 4.6). Comparing 
with the results obtained by [Cumberbatch, Lattanzi fc Silk] 
(20101, we notice that for their different SUSY neutralino 
models, the authors obtain ASTt curves that cover a simi- 
lar dynamic range than our curves. However, in the redshift 
range of 20 < z < 46 that is interesting for detection, they 
tend to predict a stronger signal than our ITeV WIMP with 
fabs = 0.01 case, to which we should compare their results 
to (because the authors already incorporate an explicit fabs 
treatment and they have values fabs ~ 0.01 in that redshift 
range; see their figure 12 and figure 13). This may be simply 
because they use a lower mass 50GeV WIMP, but they ob- 
tain the energy and number of photons and electron-positron 
pairs that result from the cascades produced by the primary 
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products of the annihilation and use these in order to com- 
pute fabs ■ Their results in the interesting redshift range are 
however similar to our fabs = 1 maximal absorption case 
with a ITeV WIMP. According to their results, all of their 
SUSY models would be in principle detectable, except for 
their model 4 (see their Table 1 for more details). 

Finally, as future work, we should mention that it is 
in principle possible to use the DM annihilation luminosity 
received by a halo (Section 5), together with the luminosity 
generated by the halo itself (Section 4.2), in order to com- 
pute the effects that the extra energy injection will have on 
the standard physical processes that take place in the bari- 
onic component of the halo (Intra Cluster Medium or ICM 
in the case of galaxy cluster-sized haloes), like for example 
its neutral hydrogen abundance, Bremsstrahlung luminosity, 
etc... Also, being able to compute the ionization fractions 
and temperatures of individual haloes at high redshifts, as 
well as their spatial clustering, may provide a way of com- 
puting the angular fluctuations of the differential brightness 
temperature {5Tb) signal for different WIMP scenarios, in 
order to have another test besides the predicted global sig- 
nal. These issues will be explored in a future paper. 
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